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Abstract 

Established results for the quark propagator in Landau gauge QCD, together with 
a detailed comparison to lattice data, are used to formulate a Poincare-covariant 
Faddeev approach to the nucleon. The resultant three-quark amplitudes describe the 
quark core of the nucleon. The nucleon's mass and its electromagnetic form factors 
are calculated as functions of the current quark mass. The corresponding results 
together with charge radii and magnetic moments are discussed in connection with 
the contributions from various ingredients in a consistent calculation of nucleon 
properties, as well as the role of the pion cloud in such an approach. 
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1 Introduction 



The structure of the nucleon has been the object of many theoretical investiga- 
tions for decades, see e.g. pp. Although many details about nucleon properties 
have emerged, some remain elusive, a fact which can easily be inferred from 
the surprises caused by recent experimental results. One prominent example is 
the discrepancy between the ratio of electromagnetic proton form factors ex- 
tracted via Rosenbluth seperation on the one hand and polarisation-transfer 
experiments on the other hand; a recent summary of the experimental sit- 
uation and possible theoretical explanations is e.g. given in Ref. [2]. While 
quark models have pioneered our understanding of the nucleon, they become 
increasingly unlikely to contribute to a more concise picture of the nucleon's 
structure: the complicated nature of baryons cannot be precisely portrayed 
by model constructions like constituent quarks or solitonic meson clouds. In 
this respect it is worth noting that efforts to bridge the gap between studies 
of strong ( "nonperturbative" ) QCD and of hadron structure are on-going and 
show encouraging results. Some examples of such progress can be found in 
the lectures compiled in Ref. [3]. Besides and together with lattice QCD and 
chiral perturbation theory functional methods have proven to be a method 
of choice when investigating the infrared properties of QCD, for some recent 
brief overview see e.g. [I]. The investigation to be reported here will build on 
such functional methods in the context of Dyson-Schwinger equations, and we 
will relate our results to corresponding ones obtained by a synthesis of lattice 
and effective field theory calculations. 

Dyson-Schwinger equations (DSEs) present a nonperturbative continuum ap- 
proach to QCD. Each Green function in QCD satisfies an integral equation, 
which involves also higher Green functions; this leads to an infinite tower of 
coupled integral equations (see e.g. Refs. [5115117115] for recent reviews of the sub- 
ject). For numerical hadron studies in QCD one relies on a truncation of the 
infinite tower of equations to a subset which is solved explicitly. Our study 
focuses on the DSE of the quark propagator and the Bethe-Salpeter equations 
(BSEs) for quark-quark and quark-diquark systems, respectively. The ampli- 
tudes obtained from the solution of the nucleon BSE are used in subsequent 
calculations of electromagnetic nucleon properties. 

Originating from the NJL model [9ff0] . the covariant quark-diquark descrip- 
tion has often been applied to the investigation of nucleon and A properties 
|11|I12|I13| . Nucleon form factors were obtained either by employing ansatze for 
the quark-diquark amplitudes [T3|T5] or by solving their BSE using ansatze 
for the ingredients of its kernel, i.e., the dressed-quark propagator, diquark 
propagator and diquark amplitudes |16|ll7lll8|ll9f20ll21ll22] . The present work 
extends these studies insofar as those ingredients are obtained consistently 
from dynamic equations at a more fundamental level of the approach. It is an 



2 



ab-initio calculation in the sense that after specifying a truncation and fixing 
the pseudoscalar meson sector to the experiment as well as interaction param- 
eters to lattice data of the quark propagator, no further observables need be 
used as an input of the calculation. 

To make this evident: the main objective of the studies to be reported here 
is to develop a QCD-based understanding of the nucleon structure. Given the 
experimental results, or more precisely, the complicated nature of the nucleon, 
the aim is as ambitious as the potential outcome is rewarding. 

This article is organized as follows: Section [2] gives a brief introduction to 
bound state equations in quantum field theory, in particular for a three-quark 
state. Section [3] introduces the quark propagator and the quark-gluon inter- 
action, i.e., the dynamical setup of the approach. The concept of diquarks 
is introduced in Sec.H] and applied to the nucleon as a two-body system of 
quark and diquark in Sec. [51 The nucleon's electromagnetic form factors in 
this framework are developed in Sec. [6) Section [7] contains our results for the 
nucleon mass and electromagnetic properties as well as a discussion; we con- 
clude in Sec. [SJ 

In order to make this paper self-contained as well as to provide all details 
necessary for reproducing our results we have added appendices containing 
solution strategies for the QCD gap equation in App. [AJ all necessary detail 
on Bethe-Salpeter amplitudes of mesons and diquarks in App. [B] as well as 
quark-diquark systems in App. £3 the construction of the electromagnetic 
current in App. [Dj and limitations on the numerical implementation coming 
from singularities of the Green functions in our approach in App. [Ej These 
also serve to avoid confusion of different notations and conventions in the 
literature. 



2 Bound-state equations in QCD 

As bound states of quarks, hadrons correspond to poles in the quark 4- or 
6-point functions and or their amputated connected parts, the scat- 
tering matrices and defined via 

G (n) = G (n) + Q(n) T (n) G H (1) 

Up to antisymmetrization, Gq is the product of n dressed quark propagators 
Si. We work in Euclidean momentum space, i.e., we use {7^,7^} = 25^ and 
y4 — anc i dropped all Dirac, color and flavor indices of each quark leg 
in the above quantities for brevity. The product in Eq. ([T|) is understood as 
a summation over all occurring indices as well as integration over all internal 
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momenta. At the pole corresponding to the bound-state mass M, bound-state 
amplitudes \I/ are introduced as the residues of the scattering matrix: 



P ,<y. 

P 2 + M 2 ' [) 

where P is the total momentum of the n quarks. The scattering kernel 

(n > 2), constructed out of /-quark irreducible components (I = 2...n), is 

related to the scattering matrix T^ 1 ' via Dyson's equation: 

T = K + KG K + KG KG K + ■■■ = K + KG T 

^ T- 1 = K- 1 - Go . [ ) 

Inserting the pole ansatz and comparing the residues leads to a bound-state 
equation and a canonical normalization condition at the bound-state pole 
P 2 = —M 2 which completely determine the amplitudes on the mass shell (see 

Fig. CD 

^ = K G ^ , (4) 
- d (K- 1 -Go) 

Eq. (j4"l) is a homogeneous integral equation and fully relativistic: it is the 
BSE in the two-body case and its quantum-field theoretical analogue for a 
three-body system. It can be solved if the quark propagator Si and the kernel 
are known. Both ingredients can in principle be obtained from the in- 
finite coupled set of Dyson- Schwinger equations (see, e.g., [8] and references 
therein). Feasible present-day numerical DSE solutions include 2- and 3-point 
functions within certain truncations, but the complexity of DSEs for 4-point 
functions has so far prevented a direct numerical approach. The construction 
of appropriate kernels is restricted by the underlying symmetries of the the- 
ory. In quantum-field theory symmetries are implemented by Ward-Takahashi 
(WTIs) and Slavnov- Taylor identities (STIs) which relate different n-point 
functions to each other. The most prominent example is the axial-vector Ward- 
Takahashi identity (AVWTI) [23], which relates the two-quark kernel to 
the kernel of the quark DSE. It is imperative to satisfy these identities in the 
truncation used in a numerical study; in this respect, the rainbow-ladder trun- 
cation has emerged as the lowest order in a symmetry-preserving truncation 
scheme pipH] . 

For baryons, the 3-body kernel in the relativistic 3-body equation can be 
written as the sum of a 3-quark irreducible contribution and the 2-body kernels 
= ® S' 1 [TTfl3f2b , | . where the subscript i identifies the spectator 
quark: 

^ {3) = ^ + E^ 3) - (6) 

i=l 
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P =-M„ 




P =-M„ 




Fig. 1. The pole ansatz for the 2- and 3-quark scattering matrices and the resulting 
bound-state equations. 

Neglecting the 3-body irreducible part yields relativistic 3-body equations of 
the Faddeev type, commonly referred to as relativistic Faddeev equations (see 

Fig. ED 

* = with % = KP ) G$ ) * = K? ) S j S k * , (7) 

i=i 

where k} denotes a cyclic permutation of {1, 2, 3}, i. e. of all Dirac, color, 
and flavor indices. The resulting relativistic Faddeev amplitude does therefore 
not contain explicit three-quark correlations. 



+ 



+ H 



Fig. 2. The relativistic Faddeev equations, Eq. (f7|). 



3 Quark propagator and quark-antiquark scattering kernel 



3.1 Elementary QCD Green functions 



Equation (JTj) is the starting point of this work. Its solution requires know- 
ledge of the quark propagator and the two-quark kernel, which are given in 
this section. The renormalized dressed-quark propagator S{p,p) contains two 
dressing functions a v and a s which correspond to a general fermion propaga- 
tor's vector and scalar Lorentz structures. Equivalently, one can write S in 
terms of the quark renormalization function A and the quark mass function 
M, which is independent of the renormalization point fi: 

S(p,fi) = (A(p 2 ,u. 2 ) (^ + M(p 2 ))) _1 = -#<7„(pV) + °s(p\n 2 ) ■ (8) 
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These dressing functions are obtained by solving the quark Dyson-Schwinger 
equation, also referred to as the QCD gap equation (see Fig. [3]) 



S-\p, fi) = Z 2 (u, 2 , A 2 ) (if + M(A 2 )) + S(p, /z, A) , (9) 



where ^(/x 2 , A 2 ) = A(p 2 = A 2 ,/i 2 ) is the quark renormalization constant and 
A is a translationally invariant ultra-violet regularization scale. The quark 
self-energy £ is defined via 

A 2 ,74 

S(p, /U ,A) = -^-Z 1F ( /U 2 ,A 2 ) l—^-irS(q,fi)D^(k,^(q,p,fi), 

(10) 

where the prefactor (JV^ — l)/(2Nc) = 4/3 stems from the color trace. £ 
contains one bare (gZ 1F ry M ) and one dressed (<?r M ) quark-gluon vertex and 
the gluon propagator D^ u with gluon momentum k = p — q which have to be 
either known a priori or determined in the course of a self-consistent solution of 
the DSEs of the quark- and gluon-propagators together with the quark-gluon 
vertex. Technical details of the solution of Eq. (j5J) are sketched in App.lAl 



-1 -1 



o 
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Fig. 3. The quark DSE in pictorial form. 



In Landau gauge, the dressed- gluon propagator characterized by a dressing 
function Z is transverse with respect to the gluon momentum k = p — q: 

D^(k,ri = ^J^Tr, Tr-.= ^-^. (ii) 

The general dressed quark-gluon vertex consists of 12 tensor structures, 

r"(p, q,fi) = f: Up 2 , q 2 , k 2 , /i 2 ) T?(k, I) , (12) 

i=l 

with the simplest possible representation of the Dirac basis elements given 

by T ?(k, i) = 7"{i, % h i th i v {i, % h i $\ and ^"{1. h h f where 1 = 

— {p + q)/2. Both gluon propagator and quark-gluon vertex satisfy their own 
DSEs. Progress on a consistent solution of this coupled system has been 
made both analytically in terms of infrared exponents of the Green func- 
tions [27,28,29], as well as numerically for general momenta in certain trunca- 
tions [30|3T] . However, the implementation of such a scheme in our problem 
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is clearly beyond the scope of the present study. Therefore we use underly- 
ing symmetry properties of QCD to arrive at a consistent truncation whose 
numerical implementation is feasible in our approach to the nucleon. 



3.2 Effective quark- gluon interaction 



The correct implementation of chiral symmetry and its dynamical break- 
ing is guaranteed if the AVWTI is satisfied, which is achieved, e.g., in the 
rainbow-ladder truncation. Using the STIs in Landau gauge: Zip = Z 2 /Z 3 
and Z g Z 3 Z 3 2 = 1, where Z 3 , Z 3 and Z g are ghost, gluon and charge renor- 
malization constants, the quark self-energy is truncated via the replacement 

g 2 Z 1F D tlv {k, //) V u (q, p, //) - Z 2 Tf 1 vf , (13) 

which leaves the effective interaction a(k 2 ) to be determined as the single 
unknown function in the approach, a is a combination of the gluon dressing 
function and a purely /c 2 -dependent dressing of the vector part of the quark- 
gluon vertex: 

a(k 2 ) = f-^Z(k 2 ^ 2 )T 1 (k 2 , f , 2 ). (14) 
4?r Z 2 Z 3 

From Z g Z 3 Z l J 2 = 1 and the renormalization-scale dependence of these quan- 
tities, given by g ~ 1/Z g , Z ~ l/Z 3 , Ti ~ Z 2 /Z 3 , it follows that a{k 2 ) is 
independent of the renormalization point. Asymptotically, it has to approach 
the perturbative behavior of QCD's running coupling, 

a{k > ' ^yiho ' (15) 

where 7 m = 12/ (llA^c — 2A^/) is the anomalous dimension of the quark prop- 
agator (in our calculation we use Nf = 4). For k 2 — > 0, an infrared analysis 
of the DSEs leads to the power laws Z(k 2 ) -> (k 2 ) 2K and T^k 2 ) -> (k 2 )- 1 / 2 ^ 
[291132] with 0.5 < k < 0.7 [27.28.33j. Accordingly, the coupling, as infered 
from Yang-Mills vertex functions becomes constant. On the other hand, the 
interaction should exhibit sufficient strength at small gluon momenta to en- 
able dynamical chiral symmetry breaking and the generation of a constituent- 
mass scale for the quark. This translates into notable non-perturbative en- 
hancements of the quark dressing functions A(p 2 ) and M(p 2 ) at infrared mo- 
menta, see, e.g., [5]. Several models for a combining the UV limit with an 
ansatz in the infrared have been employed in the past and applied to detailed 
studies of meson physics p4,35,36,37,38,39j40]. In our present study we start 
from the interaction of Ref. [39] which was fitted to light pseudoscalar meson 
data and provides an efficient description of pseudoscalar and vector mesons 
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In an attempt to combine Dyson-Schwinger and lattice approaches, several 
parametrizations for the quark-gluon interaction have been tested which are 
capable of reproducing recent lattice data at higher quark masses [54T55] for 
the quark dressing functions M and A [56"(57II58| . In the same way we tune 
the parameters of the effective interaction of Ref. [39J such that the resulting 
mass function agrees with the quenched lattice results from Ref. [51]. The 
parametrization of the interaction reads 

/f2 , C £ 7T (k 2 \ 2 k 2 /(w 2 A 2, 7T 7m (l - e- fc2 / A o) 

a (k 2 ) = -L- — e - k A ° ] + p ^ '- . 16 

-I \W Iln{ e2 -l + (l + P/A^) 2 } 

The first part of ffl6|) provides the characteristic infrared strength while the 
second term accounts for the ultraviolet behavior of f[T5"j) . Furthermore, we use 
A.qcd — 0.234 GeV and A = 1 GeV. We note that in contrast to the expected 
infrared behavior this ansatz behaves as a(k 2 )/k 2 — > const, for k 2 — > 0. This 
facilitates the numerics as the self-energy integral is of lesser divergence (for 
details, see App. [A]), but does not play an important role, since the quark DSE 
is dominated by the interaction at intermediate momenta [59]. For the light- 
quark mass, all parameters have been chosen according to the original values 
of the one-parameter model of Ref. [3S] to reproduce the phenomenological 
quark condensate and meson observables at the u/d quark mass. There and 
in subsequent meson studies, the parameter uj is used to vary the width of the 
IR part of the interaction, while the infrared strength is kept the same, which 
always leads to the correct amount of dynamical chiral symmetry breaking. In 
our present study, we fixed also the value of uj via the comparison to lattice 
data for the quark propagator. In order to match the lattice data available 
at larger quark masses (cf. Fig.HJ), both coupling strength c £ and coupling 
width uj£ must diminish with increasing quark masses. This anticipates to some 
extent a further quark-mass dependent structure in the quark-gluon vertex 
which is missing in rainbow-ladder truncation. Our fit ot data in Ref. [53] 
provides the parametrizations c £ = 0.443/(1 + 4.079 £ — 1.026 £ 2 ) and oj^ = 
0.377 + 0.295 c £ with the quark mass parameter £ = M(2.9 GeV)/AgcD- For 
the fitting procedure only we have replaced our usual value of Nf = 4 by 
Nf = 0. The renormalization scale fi = 2.9 GeV has been chosen to allow for 
an optimal fit to lattice data. In these lattice calculations the scale has been 
fixed by requiring the string tension to be a = (440 MeV) 2 . 



3.3 qq and qq bound states 

The solution of the homogeneous BSE is the on-shell Bethe-Salpeter amplitude 
(BSA) r(p, P) for a quark-antiquark bound state with relative momentum p, 
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Fig. 4. Quark dressing functions: Quenched lattice data [53] vs. the DSE result 
employing the interaction (|16p with Nf = 0. The quark masses in the legend are 
given at p 2 = (2.9 GeV) 2 . Similar results using different parameterizations of a(k 2 ) 
have been obtained in [56)57)58] . 



total momentum P, and mass M (at P 2 = — M 5 



A 



r«tf(p,P) = f-^K ai>sp (p,q,P){S(q + )T(q,P)S(-q„)}^ , (17) 

where Greek indices represent Dirac, color and flavor indices, and the quark 
and antiquark momenta are q + = q + aP and g_ = —q + (1 — a) P. The am- 
plitude's dependence on p, P can be formulated in terms of the three Lorentz 
invariants p 2 , P 2 , and p ■ P, where the latter depends on the angular variable 
z :— p ■ P I \Jp 2 P 2 . Calculations are simplified if one expands the scalar am- 
plitudes in T in Chebyshev polynomials in z (for details, see App. [B]). Due to 
Lorentz invariance, the momentum partitioning parameter a G [0, 1] can be 
chosen arbitrarily; for equal quark and antiquark masses a value of 1/2 techni- 
cally simplifies the BSA. The arguments q 2 ± of the quark propagator dressing 
functions in (fTTl) are complex for timelike P 2 ; methods to evaluate the quark 
propagator in the complex plane of Euclidean four- momentum-squared are dis- 
cussed in App. [A] The kernel K is the amputated quark-antiquark scattering 
kernel which is irreducible with respect to a pair of qq lines. The kernel con- 
sistent with the rainbow-truncated kernel of the quark self-energy is a ladder 
kernel and written as 

K a7 ,es(p,q,P) = Z\ (y) ^(f d-fUTT (vf)tu , (18) 

where the SU(3)c Gell-Mann matrices have been explicitly stated (in our case, 
the flavor structure in the BSE gives no contribution since we consider an 
equal-mass system and assume isospin symmetry). This truncation describes 
an iterated dressed-gluon exchange between the quark and the antiquark and 
respects the AVWTI and chiral symmetry: in the chiral limit it ensures that the 
ground-state pseudoscalar mesons are the massless Goldstone bosons linked 
to dynamical chiral symmetry breaking, while for finite current quark masses 
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Table 1 

Mass, strength and width parameters in the coupling Eq. (|16[) at the physical point. 
The meson masses were obtained by solving pseudoscalar and vector meson BSEs, 
the pion decay constant from Eq. (jB. 1T|) . 



m. 



0.036 0.387 0.491 138 MeV 92MeV 727 MeV 

Table 2 

Possible constituent-quark mass scales: M(0), M E defined by M§ := M 2 (p 2 
and m 

see also the 



P 



i q — |IrnyPg ing |, where p 2 ing is the location of the singular!! 

of complex conjugate poles) in a v , a s which is closest to the origin 
discussion in App.[Ej The values for the current quark mass and quark condensate 
at 2 GeV were obtained by evolving the renormalization-point independent current 
mass and condensate as determined from the perturbative tails, Eq. (|A.8|) . according 
to the 1-loop formulas employing the quenched scale Aqq D = 0.225 GeV |57|60j . 



■w 



MS 
2 GeV 



M(0) 



M E 



m„ 



MS 
2 GeV 



4.8 MeV 0.37 GeV 0.32 GeV 0.52 GeV -(255 MeV) 3 

it leads to a generalization of the Gell-Mann-Oaks-Renner relation [23] . The 
spin structure of the amplitudes for different sets of quantum numbers and 
the method for solving the BSE are presented in App.lBl 
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Fig. 5. The pion decay constant f n determined from Eq. (|B.17j) vs. quenched lattice 
results [61 1 . 



By employing the rainbow-ladder truncation of the DSEs, we link our study 
of nucleon properties to sophisticated, consistent meson studies. Our results 
for diquark amplitudes are connected directly to corresponding meson results 
in this truncation which has provided a good description of pseudoscalar and 
vector meson observables. Tabled] shows the parameters appearing in the ef- 
fective coupling at the physical light-quark mass, defined by = 138 MeV, 
together with basic meson observables. Further tensor structure in the quark- 
gluon vertex (especially a scalar part ~ l v in Eq. (TT2"]) ) seems to be necessary 
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to successfully describe other meson channels [40.62.63j. Efforts towards going 
beyond rainbow-ladder truncation have been made, e.g., in [59 62 64 6 5|66] , 
but are still numerically unaffordable in our study of the nucleon. 

By adjusting the coupling to lattice data, the current-mass (or, equivalently, 
pion-mass) dependence of meson and nucleon properties are completely deter- 
mined by the mass dependence of and uj^ in the effective coupling and there- 
fore directly relatable to quenched lattice results for the quark propagator. In 
this way, all parameters are fixed and nucleon properties are predictions in the 
model. As an example we have plotted our results for the pseudoscalar meson 
leptonic decay constant vs. corresponding lattice data in Fig. [51 In Tab. [2] we 
list numbers characterizing the quark propagator dressing functions. 



4 The diquark ansatz for the 2-quark scattering matrix 

After having specified the quark propagator, the quark- ant iquark and there- 
fore also the quark-quark kernel in terms of a rainbow-ladder truncation, the 
relativistic Faddeev equation ([7]) can be solved numerically. However, in prac- 
tice the Faddeev amplitude as a 3-particle amplitude depends on three 
four-momenta and in the case of the nucleon consists of 32 Dirac tensors, 
which renders a direct solution of Eq. ([7]) numerically expensive. 



4-1 Diquark correlations 

A possible simplification is to reduce the three-body to a two-body problem. 
The underlying assumption of the Faddeev truncation in Sec. [2] was to view 
correlations of two quarks as the dominant structure in the nucleon. This is 
motivated by the observation that colored two-quark states can appear in 
an SU(3)c anti-triplet or sextet configuration, allowing the formation of a 
color-singlet nucleon out of the anti-triplet state together with a color-triplet 
quark. In the following, we identify the dominant structures of these two-quark 
correlations in the form of diquarks, leading to a description of the nucleon as a 
bound state of quark and diquark. Further motivation of the significance of the 
diquark concept has recently come from investigations of diquark confinement 
in Coulomb-gauge QCD [BT] . 

The first step in such a construction is to rewrite the relativistic Faddeev 
equation by replacing the 2-quark kernels K\ = K\ <g> S i by 2-quark 
scattering matrices t/ 3 -* = <g> S^ 1 (i denotes the spectator quark). In 
general, they are connected to each other by the relation T" 1 = — G , 
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Fig. 6. Relativistic Faddeev equation (|20p which involves the 2-quark scattering 
matrix instead of the 2-quark kernel. 

i. e. Dyson's equation from Eq. ([3]). Then it follows from Eq. (J7J) that 

= T, t (Kr 1 - G ) * = Ti Kr 1 (tf - = (1 + T t G ) (* - , (19) 

where for the moment we omitted the superscripts in T^ 3 \ Kf^ and G^. The 
relativistic Faddeev equation then becomes a set of coupled integral equations 
for the Faddeev components (c.f. Fig. [6]) 



T^Gf^i + vl/,) = If J ^ S k (vl/, + * 



(2) 



(20) 



where the components ^ are still three-body amplitudes. For distinguishable 
quarks, {ijk} is an even permutation of {123}; for identical quarks an an- 
tisymmetrization in the total Faddeev amplitude \I> is required. Determining 

(2) (2) 

the two-body T-matrices T/ from the corresponding two-body kernels K\ 
according to Eq. (131), e.g. by inserting rainbow-ladder kernels, would make 
Eq. (120]) simply an intricate way of reexpressing the original equation (J7J). In- 
stead, we follow a different path and aim at an ansatz for the T-matrix which 
is (a) based on the assumption that it contains diquark poles at timelike values 
of the total momentum P 2 , but (b) still retains some of the characteristic fea- 
tures of Eq. ([3]). This is realized as a separable sum over diquark correlations. 
We restrict ourselves to the + and 1 + channels, i. e. to scalar and axial- vector 
diquarks for reasons explained below: 



T ay ,f3s(p, q, P) 



£r^,p)z^)(p 2 )ri;^ 



-p) 



(21) 



(pu) 



r^g are diquark amplitudes analogous to the corresponding meson amplitudes 
and the D^ v ' denote diquark propagators. The assumed poles of the T-matrix 
in the ansatz (12TI) are embedded in the diquark propagators and define the 
corresponding diquark masses: 



D{P 2 



-Mr 



P 2 + Ml 



D^(P 2 ) 



-Mr 



1 P 



P 2 + ML 



(22) 



where Tp V again denotes a transverse projector. In the same manner as de- 
scribed in Sec.|2] this leads to a homogeneous diquark BSE, T = K^G^T, 
for a diquark bound-state on the mass shell P 2 = —M 2 which resembles the 
meson BSE, Eq. ifPTj) . and was used for detailed studies of diquarks, e.g., in 
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[68,69]. The equation is shown diagrammatically in Fig. [7] and reads 

rg(p,P) = j-^r 4 K aJ MP,Q,P) {S(q + )T^( q ,P)S T ( q ^}^ , (23) 



where the momenta have been defined in the discussion of Eq. ( Tl7|) . The re- 
placements S(-q-) — > S T (q-) and K ai jp — > K ai ,ps originate from substitut- 
ing an antiquark- with a quark leg. For the sake of consistency the kernel K has 
to be identified with the rainbow-ladder kernel (fl8|) . By working out the color 
factors (see App. lB.l]) one finds that the resulting equation for V with the 



quantum numbers J p is identical to the equation for a color-singlet J~ p me- 
son except for the diquark's coupling strength which is reduced by a factor of 
2 = N c — 1. This confirms that the interaction in the color anti-triplet diquark 
channel is strong and attractive. The same analysis also shows that the inter- 
action is strong and repulsive in the color sextet channel [TOfTI] . Comparison 
with meson phenomenology therefore suggests that the lowest-mass diquarks 
are the scalar diquarks (the partners of the pseudoscalar mesons), followed 
by axial-vector, pseudoscalar and vector diquarks. This was also observed in 
Bethe-Salpeter [68] and lattice [72] investigations and justifies the restriction 
to scalar and axial-vector diquarks for describing light baryons composed of 
quark and diquark. 



P- q- 




Fig. 7. The Diquark BSE (l23l). 

The ansatz (I2T]) is useful only if the full 2-quark T-matrix, as obtained from 
the quark-quark kernel via Dyson's equation ([3]), indeed contains timelike 
scalar and axial-vector diquark poles. Since asymptotic colored diquark states 
correspond to timelike poles in the diquark propagators, one could suspect 
a violation of diquark confinement in this case. However, the absence of a 
Lehmann representation of a certain propagator is a sufficient but not neces- 
sary criterion for confinement of the corresponding state due the associated 
violation of reflection positivity [6l73|74j . In this way, two-point correlations 
of colored fields may contain real timelike poles in momentum space without 
contradicting confinement, a statement which is also true for the quark prop- 
agator [33ll75f76j . Free-particle quark and diquark propagators can also yield 
quantitatively meaningful results for hadronic observables (see, e.g., [19]). 
The existence of a solution of (1231) implies that a rainbow-ladder truncation 
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Table 3 

Scalar-axialvector diquark mass splitting in the chiral limit. The BSE value is com- 
pared to several lattice-QCD results. The units are GeV. 

BSE Ref. [78J Ref. [72] Ref. p] Ref. [80] 
0.21 0.10(5) 0.14(1) 0.29(4) 0.36(7) 

induces such diquark poles, and for our particular approach and methods this 
pole structure is in fact a necessary prerequisite. On the other hand, the intro- 
duction of interaction terms beyond rainbow-ladder truncation in the skeleton 
expansion of the quark-quark kernel in the diquark BSE removes diquark 
states from the physical mass spectrum because of large repulsive corrections 
[24l25f64|77] . Nevertheless, also kernels that do not produce diquark bound 
states support a physical interpretation of the poles in the diquark propaga- 
tors either in terms of mass scales or inverse correlation length scales of the 
diquark correlations inside a baryon. 
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Fig. 8. Scalar and axial-vector diquark masses together with the quark constituent 
mass scale M(0) vs. squared pion mass. The vertical line marks the physical point. 

Diquark masses are gauge- dependent, but gauge-independent mass differences 
can be determined from lattice calculations. Several such investigations have 
been performed over the last years in different approaches and with various 
fermion actions [72|78f79f80f81j . While they exhibit quite different quantitative 
results, the common qualitative feature is that the mass splitting between 
scalar and axial- vector diquark in the chiral limit is of the size of hundred to 
several hundred MeV (see Tab. [3]) and decreases with increasing current quark 
mass; Fig. [8] shows scalar and axial- vector masses together with a measure of 
the quark constituent mass as result of the DSE-BSE approach employed in 
this work. 
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4-2 Offshell and asymptotic behavior of the T-matrix 



By solving the scalar and axial- vector diquark BSEs, Eq. (|23|) . the quark-quark 
scattering matrix is known at the diquark mass poles, i. e., at P 2 = — M 2 C and 
P 2 = —M 2 V . The description of baryons as composites of quark and diquark 
also requires the knowledge of the T-matrix for general diquark momenta. The 
ansatz (1211) dictates its off-shell behavior via both off-shell diquark amplitudes 
and diquark propagators. Inverting this ansatz and inserting Dyson's equation 
([3]) determines the inverse scalar and axial-vector diquark propagators for all 
values of P 2 (the superscripts in T^ 2 \ and Gq are omitted): 



Eq. ([21]) necessitates knowledge of the diquark amplitudes' off-shell behavior. 
This cannot be determined without further information; a possible ansatz is 
constructed in App. IB.3l 

From Eq. ([3]) one infers that T is asymptotically dominated by the gluon ladder 
exchange which is constant in P 2 : the ladder diagram is the lowest term in 
a skeleton expansion of K and independent of P 2 (all higher contributions 
drop off for P 2 — > oo) and, as a product of two quark propagators, G ~ 
P~ 2 for P 2 — > oo. Eq. (I2ip . as a sum of two separable terms, is probably no 
accurate representation of the ladder kernel in the ultraviolet, but the correct 
asymptotic power behavior T DT — > const, is guaranteed via Eq. (1241) . If the 
diquark amplitudes asymptotically behave as P x , then D^ 1 — ► P 2X such that 
T becomes constant in the ultraviolet. Eq. (1241) is explicitly worked out in 

App.E! 



5 The quark-diquark BSE 

All elements of the relativistic Faddeev equation (120]) . i.e., the dressed-quark 
propagator and the quark-quark scattering matrix, have been specified in 
the previous sections. The separability of the diquark ansatz (1211) for the T- 
matrix allows for the reduction of the original three-body equation for the 
nucleon to a two-body Bethe-Salpeter equation for a quark-diquark bound 
state. Up to symmetrization, quark-diquark amplitudes are introduced by re- 
moving a diquark amplitude and propagator from the Faddeev component 




(24) 



# a/37 of Eq.®: 



*afh(p,Pr,P) = J2 T h(PnPd)D ab ( Pd ) {&(p, P) u(P) } q . (25) 



a,b 
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a, b, c are diquark indices: T 5 denotes the scalar and r a=1 " 4 the axial- vector 
diquark amplitude; the diquark propagator D ab is either scalar (a, b = 5) or 
axial- vector (a,b = 1 ... 4). The momenta p and p r are the relative momenta 
between quark and diquark and within the diquark, p d and P are total diquark 
and nucleon momenta. The spinors u a (P) are solutions of the Dirac equation 
in Euclidean space. The remaining matrix-valued Bethe-Salpeter amplitudes 
$a*(p, P) contain one (scalar or axial- vector) diquark and two fermion legs 
and feature a decomposition constructed from the same Dirac basis elements 
as used in the meson and diquark case (see App.[U]). The spin- and isospin- 
1/2 nucleon is therefore a sum of scalar and axial- vector diquark correlations, 
and the quark-diquark amplitude describes the relative momentum correlation 
between quark and diquark. Inserting the ansatz (I2T!) for the T- matrix together 
with (125]) into the relativistic Faddeev equation (1201) yields a quark-diquark 
Bethe-Salpeter equation 



Kb(p,p) 



d A k 



{K ab (p, k, P) S(k q ) <l> c (k, P)} ap D bc (k d ) 



' J (2tt) 4 
where the quark-diquark kernel is given by 

K%[p, k, P) = {T\k r , k d ) S T (q) T a (p r , -p d )} 



a3 



(26) 



(27) 



which couples scalar and axial- vector diquark amplitudes (i.e., a, b = 1 ... 5). 
It describes an iterated exchange of roles between the spectator quark and 
the quarks which constitute the diquark; this quark exchange generates the 
attractive interaction that binds quarks and diquarks to a nucleon. 





Fig. 9. The quark-diquark BSE (126ft in pictorial form. 



The momenta in Eqs. ( 12611271) are (cf. Fig.[9]): 



Pq=P + P , 
kg = k + 7] P 



Pd= -p+ (l-rj)P , 
k d = -k + (1 - 77) P , 



q=Pd~ k q , 

k r = (l-a)p q -aq, (2S 
p r = (1 - a) k q - aq . 



Again, the momentum partitioning parameters a, 77 G [0, 1] for diquark and 
quark-diquark amplitudes are arbitrary since in a covariant description there 
is no unique definition of a relative momentum. Translation invariance implies 
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that, e.g., for each BSE solution $(p, P, rf) a family of solutions of the form 
$(p + (?/ — r))P, P,r)') exists |20j. We set a = 1/2 but keep r\ as a variable since 
it can be used to ease the constraints caused by singularity structures in the 
complex plane (cf. App.[E|). By working out the color and flavor factors of the 
quark-diquark amplitudes, given in Eq. fIC.ll) . and of the diquark amplitudes 
(IB.6tlB.7|) . the BSE kernel picks up a color- flavor factor 



where the first row (column) represents the scalar part of the kernel and the 
second row (column) the axial- vector part. 

Eq. (126]) is solved using similar methods to the meson- and diquark-BSE cases, 
in particular: (a) reduction to a system of homogeneous one-dimensional inte- 
gral equations by a Chebyshev expansion of the coefficient functions in their 
angular variable, and (b) introducing an artificial eigenvalue which becomes 1 
at the bound-state mass. These procedures are explained in detail e. g. in (82]. 
Solving the quark-diquark BSE yields the nucleon mass and quark-diquark 
amplitudes on the nucleon's mass shell. Since the calculation can in principle 
be performed at arbitrary current masses and the corresponding pion mass is 
easily obtained, one can determine the nucleon mass as a function of the pion 
mass and directly compare it with lattice investigations and chiral extrapola- 
tions. Such results are presented in Seed 



6 Electromagnetic current 

The construction of an electromagnetic current operator in the framework of 
Bethe-Salpeter equations was first treated by Mandelstam [83J. In our ap- 
proach, incoming and outgoing nucleon states are described by quark-diquark 
amplitudes. In terms of the interaction with an external current, the baryon is 
resolved into its constituents, quark and diquark and the interaction between 
them, to each of which the current can couple [17]. A systematic procedure 
for the construction of the nucleon-photon vertex based on electromagnetic 
current conservation is the "gauging of equations" prescription [84, 85, 86]. In 
this context, "gauging", formally denoted by T —>■ T M , is a derivative: it is 
linear and satisfies Leibniz' rule. For a general n-body bound-state amplitude 
\l/ the current matrix is defined as the residue of the gauged scattering matrix 
at the bound-state pole with mass M: 




(29) 



P 2 ^-(M 2 +Q 2 /4) 



(30) 



T 



P 2 + M 2 



(P 2 + M 2 + Q 2 /4) 2 
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where K is the interaction kernel, T is the scattering matrix and Go the 
product of n propagators; P is the total momentum, 

P = (P t + P f )/2 (31) 

is the Breit momentum (i.e., the average of incoming and outgoing total mo- 
menta), and Q = Pf — Pi the photon momentum. From = —T (T _1 ) /i T, 
T- 1 = K- 1 - G and [K^f = —K~ l K^K~ l one finds 

T M (m2+q2/4) > T (Gff + G K^G ) T (32) 

and therefore 

= f (erg + G A^G ) * . (33) 

In our present context, $ has to be identified with the quark-diquark ampli- 
tude $, T is the quark-diquark scattering matrix, Go = S D the product of 
a dressed quark and diquark propagator, and K = T ST the quark-diquark 
kernel describing the quark exchange. The quark-photon and diquark-photon 
vertices are defined as the gauged inverse propagators: := — (S' _1 ) M and 
:= — (D~ i y ' . The gauged diquark amplitudes (or seagulls) T M =: M M de- 
scribe the photon coupling to the diquark amplitudes. The ingredients of the 
current matrix are therefore written as 

G^ = {SDf = ST^SD + SDT^D , (34) 

k*= (rsfy = M f *sr + rsr%sr + rsM fl . (35) 

These diagrams are worked out in detail in App.[D]and depicted in Fig. ID.ll 
The current matrix elements related to the gauged quark-diquark propagator 
Gq constitute the impulse- approximation diagrams where the photon couples 
to either quark or diquark line. The premise of current conservation also entails 
the necessity of including the diagrams containing the gauged kernel [17] . The 
electromagnetic current of a baryon in a quark-diquark framework is therefore 
completely specified by identifying the quark-photon vertex, the scalar and 
axial- vector diquark-photon vertices and an ansatz for the seagull terms. These 
quantities are constrained by Ward-Takahashi identities and thereby related 
to quark and diquark propagators and diquark amplitudes which have already 
been determined previously. This is demonstrated in Apps. HXTl and HX2l 

In order to establish a link between the electromagnetic current (133]) and 
electromagnetic form factors, one recalls that the most general form of the 
nucleon-photon current is obtained by sandwiching the fermion-photon vertex 
between spinors u(Pf, Sf), u(Pi, Si) which are solutions of the Dirac equation: 
A+(Pi)u(Pi, Si) = u(Pi,Si), u(Pf, Sf)A+(Pf) = u(Pf,Sf), where Pj and P f with 
Pf — P? — —M 2 are initial and final nucleon momenta, Si, Sf = ± are the 
spin labels, and e.g. A+(Pj) = £ Si u(P u Si)u(P u s<) = {1 + Pi/(iM)}/2 is the 
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positive-energy projector of the incoming nucleon. The matrix-valued current 
of Eq. (|33|) is obtained by removing these spinors via taking the spin sums, i.e. 
contracting the vertex with the projectors themselves: 



r(Q 2 ) = A + (P f )(F 1 zY - F 2 ia^Q v )A + (P i ) 

= A + (P f )((F 1 + 2MF 2 )ir-2F 2 P>*)A + (P i ) , '" >b> 



with <j iw = — i[y^, l v }/2. The second row was obtained by using the Gordon 
identity 

A+(P/) Kt^ ~ plM+ lMl ) A+{Pi) = • (37) 



The two dressing functions Fx and F 2 are the Dirac and Pauli form factors 
and depend only on Q 2 since on the nucleon's mass shell P 2 = —(M 2 + Q 2 /4) 
and P ■ Q = 0. F\ is dimensionless and F 2 has the dimension of an inverse 
mass; for Q 2 = they reduce to the nucleon's charge and anomalous magnetic 
moment. Note that we did not introduce the Pauli form factor in terms of 
nuclear magnetons, F 2 = F 2 /(2M) because the nucleon mass M employed 
here is the result of the quark-diquark BSE and usually not identical to the 
physical mass M exp = 0.94 GeV. Furthermore, M will increase when going to 
larger current masses and it is advantageous to use magnetons of a fixed mass 
scale for comparison with lattice calculations. 

The requirement Fi(Q 2 = 0) = 1 for the proton is automatically satisfied if 
the quark-diquark amplitudes are canonically normalized via Eq. flC.61) |17j . 
The Sachs electric and magnetic form factors, defined as 

Ge = Fx — 2MtF 2 , 2MG M = F x + 2MF 2 , (38) 

with r = Q 2 /(AM 2 ), can be extracted from the current fl36l) via 

G e = ^Ltt {rP»} , 2MG M = ^ If Tr {J^} ■ (39) 

Given a form factor F(Q 2 ), the corresponding electromagnetic radius r> is 
defined as the slope at zero momentum transfer via the Taylor expansion 



(40) 



F(0) is the corresponding electric charge or magnetic moment which in the 
case of the neutron's electric form factor, G E (0) = 0, is omitted in Eq. (liDjl . 
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7 Numerical results and discussion 

To facilitate the discussion, we shortly recapitulate the steps employed in 
this work. Starting from the general relativistic three-body equation for the 
baryon's amplitude PJ, we neglected the three-body irreducible kernel, assum- 
ing that correlations of two quarks are dominant in the nucleon, and arrived 
at a relativistic Faddeev equation ([7]), which we converted to (l2"Ui) such that 
it involves the two-quark scattering matrix T. Instead of determining T from 
Dyson's equation ([3]), we employed a separable ansatz fl2Tl) involving diquark 
amplitudes and a diquark propagator which exhibits timelike poles at scalar 
and axial-vector diquark masses. With this simplification Eq. (120]) could be 
reformulated in terms of a quark-diquark Bethe-Salpeter equation f)2"o]) for the 
baryon amplitude. The ingredients of this equation are: 

• The dressed-quark propagator determined from its Dyson- Schwinger equa- 
tion via a rainbow truncation. We employed a parametrization for the 
quark-gluon interaction (IToT) which was adjusted to basic meson observables 
and to quenched lattice results for larger current-quark masses. 

• A diquark propagator that is obtained from Eq. (1241) . which immediately 
follows from the ansatz for the T-matrix and needs as an input the 

• diquark amplitudes obtained by solving a diquark BSE (TTTT) at the poles 
in the T-matrix. For offshell diquark momenta, an ansatz is employed for 
these amplitudes (IB.20HB.23I) . 

The quark-diquark BSE yields quark-diquark amplitudes which are needed 
to calculate the electromagnetic current diagrams ()33H35p . They depend on 
the quark-photon vertex, the diquark-photon vertices and the seagull ver- 
tices. For the quark-photon interaction we employed the Ball-Chiu vertex 
(i.e., Eq. (1D.20I) without the transverse part). The diquark-photon vertices 
are given by (ID. 271) and the seagulls by Eqs. (ID.34HD.l35l) . All of the subse- 
quent results were obtained by retaining only the leading (zeroth) Chebyshev 
moments of the diquark amplitudes. The necessity of this simplification is 
argued in App.[E] Concerning the off-shell dependence of the diquark ampli- 
tudes, we use n = 1 in Eq. (IB. 241) . i.e., a moderate -v/P^-like suppression of 
the non-leading amplitudes at perturbative momenta. Occasionally we will 
also compare with the choice n = 2. 



7.1 Nucleon mass and pion- cloud effects 

Table |4] shows the pion decay constant and meson, diquark and nucleon masses 
at the physical point as obtained from Eqs. (1B.17|) . (ITT)) . (123)) and (126)) . The 
results for m p and the diquark masses are comparable to the values obtained 
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Table 4 

Results for pion decay constant and meson, diquark and nucleon masses (in GeV) 
at the physical point, defined by m n = 138 MeV. "Pull" denotes the complete 
calculation and "dom." the case where only the dominant diquark amplitudes ~ 
7 5 C, 7^C are retained. The deviations compared to Table [1] are due to the omission 
of the angular dependence in meson and diquark amplitudes. 





u 


m p 


m sc 




M N 


full 

dom. 


0.101 
0.072 


0.736 
0.773 


0.676 
0.644 


0.889 
0.894 


0.931 
0.968 


exp. 


0.092 


0.770 






0.94 



in Refs. [68] and [39] for uo = 0.5 GeV (in the present work, c<jg = 0.49 at the 
physical point, cf. Tabled]). At first sight, the result M N = 0.93 GeV for the nu- 
cleon mass seems quite remarkable since after fixing the coupling strength (1161) 
no further parameters (in particular, no observables related to the nucleon) 
have been used as an input of the calculation. Matching the experimental nu- 
cleon mass therefore implies that the corrections to the truncations made in 
our calculation cancel. In particular, one would expect corrections from drop- 
ping the quark-diquark assumption (i.e., reverting Eq. ( 1261) to a relativistic 
Faddeev equation with a ladder truncation in the quark-quark channel), cor- 
rections from going beyond rainbow-ladder truncation (towards a full 2-quark 
kernel and quark-gluon vertex), and also corrections from irreducible 3-quark 
contributions. An investigation of the large quark- mass behavior (cf. Fig. fTOl) 
indicates that the cancellation of these effects at the physical point is acci- 
dental. Still, it shows that the quark-diquark picture accounts for more than 
~ 90% of the values of the nucleon mass as obtained from lattice investigations 
at different current quark masses. 

From a phenomenological point of view, an interesting contribution comes 
from the so-called pion cloud of the nucleon. In terms of quarks and gluons, 
it corresponds to the long-range part of qq correlations interacting with the 
nucleon. Established in the cloudy bag model [92|93|94] , where the pion field 
is coupled to a constituent-quark bag [95], it is considered to be an important 
component of the nucleon's structure at low energies and small quark masses. 
Pion effects can be treated in a systematic way by chiral effective field theory, 
which is an approach to describe low-energy QCD with effective pion, nucleon, 
and A degrees of freedom [961197] . In combination with lattice simulations, it 
has proven to be an efficient tool for describing masses and electromagnetic 
properties of hadrons [98H99II100] . In this framework the nucleon mass near the 
chiral limit is obtained from the chiral expansion [101] 

M N (ml) = oo(A) + a 2 (A)ml + a 4 (A)m£ + ■ ■ • + £(m£ , A) , (41) 

where a, are the (bare and a priori unknown) parameters appearing in the 
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Fig. 10. The DSE/BSE result for the nucleon mass (dashed line: dominant diquark 
amplitudes only, solid line: full calculation) versus a compilation of contemporary 
lattice data, both quenched (open circles [57]) and unquenched (filled squares [55] . 
filled diamonds [89]). Dotted [SO] and dash-dotted lines [91J with respective error 
band represent a chiral extrapolation and an interpolation between the physical 
point and selected lattice data, respectively. 

effective Lagrangian that correspond to the "nucleon core", and E denotes 
the sum of all meson-loop contributions, i. e., to 1-loop order the sum of Ntt, 
NAtt etc. self-interactions (cf. Fig. fTTl) . A is a generic regularization parameter 
that appears because effective field theory treats baryons and pions as point 
particles. This is appropriate for long-distance physics but leads to divergences 
in the loop integrals since the composite substructure is not accurately imple- 
mented [102J. A non-pointlike nucleon-pion interaction is physically equivalent 
to using a momentum cutoff corresponding to the baryon size (A > 0.2 GeV 
~ 1 fm _1 ) which picks out the long-distance or low-energy part of the self- 
energy integrals. Then, by expanding in m%, £ is split into cutoff- independent 
non-analytic terms, Eva^ ~ Tn%, Sjvatt ~ rri^lnm^, and cutoff-dependent 
terms that are even in ra£. Since the A-dependence of the latter is of the same 
type as in the bare coefficients ai, both can be combined to renormalized 
coefficients such that the nucleon mass finally reads [lOlj 



where usually the experimental numbers for the pion decay constant (f n = 92.4 
MeV), axial coupling (g^ = 1.26) and nucleon-delta mass splitting (A = 
0.292 GeV) are inserted. The final physical result is expressed in terms of 



M N {ml) = M N (0) + c 2 ml + c 4 m* + . . . 




£ In m n + . . . , 



(42) 
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renormalized low-energy constants Mjv(O), C2 and C4 which are determined 
by a fit to lattice data. Through this renormalization procedure the regulator 
dependence is removed, at least if one would work to all orders of the chiral 
expansion. To finite order a small residual cutoff dependence is left, therefore 
a convenient choice of the regulator can be utilized to improve the convergence 
of the chiral series |102j . 



Fig. 11. Expansion of the nucleon propagator in chiral perturbation theory. Solid, 
dashed and thick solid lines correspond to nucleon, pseudoscalar meson and A 
degrees of freedom. 

In order to study the (cutoff-dependent) separation of the nucleon mass into 
a "core" and a "meson cloud" contribution one has to examine Eq. (j4ip . since 
in ( 142]) the meson cloud effects in the chiral limit, E(0, A), have already been 
absorbed into Mjv(O). For a dipole regulator with A ~ 0.8 GeV, the typical 
value is of the size S(0,A) ~ -250 MeV [90ll98|ll01| . therefore leading to an 
expected core mass of ao(A) ~ 1.2 GeV. Different shapes of the regulator can 
vary this result considerably; a dipole form mimics the physical shape of the 
meson-baryon vertex, e.g., in terms of the observed axial form factor of the 
nucleon [Jj. Similar numbers have been reported from a Dyson-Schwinger so- 
lution of the nucleon propagator dressed by pions (~ —200 MeV) [T8"lll03| . in 
studies of the cloudy bag model (—300 to —400 MeV) [104J, and in a perturba- 
tive quark-diquark study with pointlike pion exchange (—150 to —300 MeV) 
[104J. This confirms that nucleon-pion loops are attractive and the binding 
energy has the effect of lowering the nucleon's mass. In agreement with the 
presumption that pion-loop effects should be suppressed in the heavy-quark 
regime, S(m^., A) is seen to fall off with increasing pion mass |98j . 




From Fig. [Til] it is clear that pion-cloud effects cannot be the only source of 
disagreement between the DSE/BSE curves and the lattice results, since we 
underestimate the data, which for large should eventually reveal the nu- 
cleon's quark core, by ~ 10%. On the other hand, a remarkable feature is 
that the slopes dM^jdm^ of both approaches are identical above a certain 
value of m^. This can be traced back to the coupling strength of Eq. ( Tl6]) 
which was fixed to (quenched) lattice results of the quark mass function. To 
evolve from the " quark-diquark core" to the full quark core of the nucleon, we 
consequently have to search for repulsive contributions that are constant in 
the quark mass. This is even more imperative when considering irreducible 3- 
quark correlations which were neglected through the Faddeev truncation (]7]): 
the leading diagram in this case is the 3-gluon vertex coupling to all three 
quark lines. It is expected to provide a further attractive contribution to the 
nucleon structure which, at least at the perturbative level, is independent of 
the current- quark mass. 
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Fig. 12. The DSE/BSE result for the p mass versus CP-PACS data [IDS] , together 
with the chiral extrapolation of Ref . [1U6] . 

It is enlightning that the same qualitative mismatch between BSE/Faddeev 
and lattice approach is already visible in the vector meson case (Fig. [T2|) where 
the rainbow-ladder truncation is the only effective restriction. The conclusion 
is that dressed gluon exchange alone binds too strongly, and that a quark-gluon 
vertex and quark- (anti) quark kernel beyond rainbow-ladder truncation should 
provide the necessary amount of repulsion in the meson and diquark channels. 
As mentioned in Sec. HI such a generalization will probably remove the timelike 
diquark poles in the quark-quark T-matrix but could still maintain a measure 
for the diquark's mass scales which enter the nucleon mass. In this respect 
we note that in Ref. [lU3j the full Dyson- Schwinger solution for the nucleon 
propagator, dressed with pions, again reduced the nucleon mass shift E(0, A) 
by increasing the width of the nucleon-pion vertex above ~ 1 GeV (loosely 
speaking: by shifting more and more short-range "pionic" contributions from 
the nucleon core into the pion cloud) . This might be a sign for the short-range 
qq correlations producing the repulsive effect that is missing in rainbow-ladder 
truncation. 

Interestingly, while there seems to be virtually no difference between quenched 
and dynamical lattice data for the nucleon mass at large m^, unquenched 
lattice calculations exhibit sizeable effects at the propagator level, e. g. for 
the quark mass function [S3]. A fit to these dynamical data probably would 
have diminished the slope dM N /dm% considerably but exceeded the parame- 
ter bounds within which the ansatz (|T6|) could be safely applied. However, the 
issue of quenching/unquenching at the rainbow-ladder level is elusive, since un- 
quenching effects in the gluon propagator alone (i. e., via solving the full gluon 
DSE) only affect the interaction a(k 2 ) in the quark DSE in terms of a slight 
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Table 5 

Proton, neutron and isoscalar magnetic moments in units of physical magnetons, i.e. 
G P M n {Q 2 = 0) = v p , n /(2M e * p ), and according to Eq. ®: fx p , n = M e * p /Mf? lc + K p , n , 
with k s = K p + K n . The electric and magnetic radii of proton and neutron are 
given in units of fm; r% denotes yj— (r^.) 2 . "Param." denotes the use of analytic 
parametrizations for the diquark-photon vertices. 





Hp 




K s 


r p E 


r" 




' M 


full 


2.52 


-1.55 


-0.04 


0.67 


0.13 


0.58 


0.57 


dom. 


2.70 


-1.68 


0.05 


0.69 


0.08 


0.60 


0.60 


param. 


3.91 


-2.01 


0.89 


0.63 


0.26 


0.45 


0.49 


exp. 


2.79 


-1.91 


-0.12 


0.87 


0.34 


0.86 


0.88 



change of its intermediate- momentum shape [107] . More drastic consequences 
for hadronic observables arise from "unquenching" the qq kernel and quark- 
gluon vertex and have been studied by implementation of an additional effec- 
tive pion exchange [59]. Finding an appropriate symmetry-preserving trunca- 
tion beyond rainbow-ladder might not only be the way to provide sufficient 
repulsive strength in the two-body channel, but possibly also the key direction 
for implementing the pion cloud in the Faddeev approach. 

7.2 Magnetic moments and electromagnetic radii 

The first and last rows of Table[5] show results for the magnetic moments and 
charge and magnetic radii together with the experimental values. The calcu- 
lated magnetic moments of proton and neutron are somewhat smaller than 
seen in experiment. All radii are sizeably underestimated, with the magnetic 
radii being yet smaller than the electric ones. On physical grounds this is ex- 
actly what is to be expected from a calculation missing contributions from the 
pion cloud. The chiral expansion for the magnetic moments in chiral pertur- 
bation theory when taking into account only the iV7r-ioop reads [108] 

^ n {ml) = //'"(0) T + • • • , (43) 

In the same manner as Eq. (1421) for the nucleon mass, Eq. (1431) does not reflect 
the cut-off dependent separation into quark core and pion contributions: the 
pion cloud gives a non- vanishing contribution in the chiral limit which has been 
absorbed into fi p,n (0). Such a decomposition has been explicitly performed in 
Ref. [109] using a dipole regulator and a cutoff A = 0.8 GeV. The result is 
displayed in Fig. [13] and exhibits a slowly decreasing magnetic moment of the 
proton's quark core, a finite and positive contribution from the pion loop in 
the chiral limit (and also at the physical point), and a finite pion contribution 
at large quark masses. The decomposition for the neutron looks similar. The 
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DSE/BSE results for the anomalous magnetic moments of proton and neutron 
are shown in the left panel of Fig.[TH The qualitatively similar behavior of both 
findings strongly suggests the long-range pion cloud to be the main missing 
contribution to the present approach. 

The right panel of Fig.[T3] displays the isoscalar combination k s = k v + K n . 
From (H3j) it is clear that the leading non-analytic parts ~ m n cancel in k s , 
i.e., up to leading order in the chiral expansion the isoscalar magnetic moment 
is not quark- mass dependent. A sizeable mass dependence of k s for small 
would therefore definitely indicate the absence of other effects apart from 
pionic corrections. The plot shows that this is not the case for small masses 
but appears above m n ~ 0.2 GeV. 

A chiral expansion similar to (1431) exists for the electromagnetic radii of proton 
and neutron [97J. In contrast to the magnetic moments the charge radii diverge 
in the chiral limit since a massless pion can propagate over infinite distances. 
The results shown in Fig. [15] provide further support for a interpretation in 
terms of a missing pion cloud as the DSE/BSE results are found to constitute 
a flat plateau at the level of the lattice data which are available at large pion 
masses. 
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Fig. 13. Meson-loop contributions to the proton's magnetic moment in Heavy 
Baryon Chiral Perturbation Theory |109| with lattice data from [87]. Picture taken 
from Ref. [109] . 



In Table [5] and Figs.fTUl [121 and [T6l we also compare with the results obtained 
by taking only the dominant scalar and axial-vector diquark amplitudes into 
account, i.e., retaining only an s wave contribution in the diquark rest frame. 
To maintain consistency throughout the whole calculation this simplification 
must be made already at the level of the diquark BSE. In this case the general 
diquark amplitudes (1B.6HB.7I) and correspondingly (1B.4I) for the meson case 
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Fig. 14. Anomalous magnetic moments of proton and neutron (left) and isoscalar 
(right). The DSE/BSE results (solid curves) are compared to quenched QCDSF 
lattice data and their chiral extrapolations (dashed curves) |100J. The dots mark 
the physical values. 
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Fig. 15. Isovector radii corresponding to the Dirac and Pauli form factors 

T?v _ rpP 
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circles [111| . filled triangles [110] ) lattice data, 
reduce to 

r(g,P) = /rV,PWc, 
r»(q,P) = fr( q \p 2 )trc. 



(44) 
(45) 



All of the previous quark-diquark studies (e.g., [T4"|15l)16Jll7Jll8l)19f20f21f22] ) 
were carried out under this truncation. Since it changes most of the consid- 
ered observables by < 10% compared to the full result (see also Ref. [68]), it 
can safely be viewed as a reasonable approximation to the overall problem, 
although it may obscure the interpretation of Mn and m p when comparing to 
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lattice data. Exceptions are the pion decay constant which is reduced by ~ 30% 
and the structure of the electric form factors in Fig. [16] around Q 2 ~ 1 GeV 2 . 



Another observation concerns the use of ansatze for the diquark-photon ver- 
tices. For instance, the general axial-vector diquark-photon vertex can be writ- 
ten in the form of Eq. (1D.22j) . of which f^ v ~ dq and f^ v dq are the dominant 
transverse components and correspond to the axial-vector diquark magnetic 
moment 

Vd q = A av ~ dq (-™1, -m 2 av , 0) = -/ 3 av ~ dq (-™L 0) . (46) 



Therefore, a possibility to circumvent the intricate vertex (1D.27|) is to use 



version (1D.22|) and neglect all the transverse terms except the one involving 
the magnetic moment, /i^ {d^Q a — S^Q^J . may be generalized to include 
an appropriate Q 2 dependence, and if only the on-shell form of the vertex is 
of interest, one can additionally contract the vertex with transverse projectors 
from the external diquark lines [21J. Reasonable values for n^q are ~ 2 — 3; the 
full quark-loop vertex (1D.27|) yields \i^ q = 2.7. We apply the same procedure 
to the scalar diquark-photon vertex where we use 

qbm and set /r~ dq to zero. 

and the scalar-axialvector transition vertex which may be expressed by 

T^(p, q) = ie^P x pPq X — , (47) 

M N 

with K sa = 2.3 extracted from the general form (ID. 271) . Table [5] shows that 
these parametrizations lead to inflated und unbalanced magnetic moments of 
proton and neutron since k s is now large and positive. The reason is that 
the on-shell values of the transversal dressing functions are not enough to 
describe the full vertex: at off-shell momenta (especially large spacelike mo- 
menta) there is a considerable amount of attenuation which is not captured by 
the above forms. The main part of the discrepancy is produced by the axial- 
vector diquark-photon vertex whose contribution to the magnetic moments is 
blown up by an order of magnitude due to the missing off-shell structure. If 
parametrizations of the above types are employed and desired to reproduce 
the magnetic moments, it is therefore mandatory to include a functional de- 
pendence on both the photon momentum Q 2 and the average in- and outgoing 
momentum P 2 . 



7. 3 Electromagnetic form factors 



Fig. [[6] shows the proton and neutron's electromagnetic Sachs form factors at 
the physical point compared to experimental results (upper part), and split 
into the most dominant contributions from which the relative strengths can be 
read off (lower part). In the case of the neutron form factor one can see that 
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Fig. 16. For the experimental points the data sets as selected in Ref. [112] are used 
[Data compiled by P. Grabmayr, Univ. Tubingen] . 

the contributing diagrams are not small by themselves: it is current conserva- 
tion that ensures the cancelation at Q 2 = such that the neutron charge is 
exactly zero. The unit charge of the proton is the result of the quark-diquark 
normalization condition (1(160 . The strongest contribution is usually the di- 
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rect coupling of the photon to the quark line. However, the shift of spectator 
quark into diquark degrees of freedom necessitates the remaining diagrams of 
Fig. ID.ll For the electric form factors it is the photon coupling to the scalar di- 
quark which is important; for the magnetic form factors the scalar-axialvector 
transition, which is a purely transverse term and does not contribute to cur- 
rent conservation, provides a small contribution. The exchange and seagull 
contributions tend to be small but are necessary to guarantee a conserved 
current. 

At larger momentum transfers Q 2 > 2 GeV 2 , pion cloud effects should vanish 
since the structure of the nucleon at small distances (< 0.2 fm) is probed. 
The large-Q 2 behavior of the proton's form factor ratio \x v G v E jG v M is there- 
fore a genuine rendition of the nucleon's quark core. Since quark and diquark 
propagators obtained from the quark DSE or via the T-matrix ansatz nec- 
essarily contain singularities, our applicable Q 2 range, depending on model 
details, is limited to below 1 — 2 GeV 2 , i.e., the region where the pion cloud 
is still effective (for a detailed discussion, see App.[E|). Access to large Q 2 
and comparison to the true "quark core" as seen in experiment can therefore 
only be established by finding appropriate methods to evaluate these quanti- 
ties beyond their dominant singularities and to include the respective residue 
contributions in the form factor integrals. 

In Fig. [T7] for the proton's form factor ratio n P G p E /G p M , the calculation involv- 
ing the dominant diquark amplitudes points towards an early zero crossing at 
Q 2 ~ 2 GeV 2 . The curve rises again when including the full diquark substruc- 
ture, and the common non- vanishing slope at Q 2 = is caused by the deviation 
between electric and magnetic radii (Table[5]), since the Taylor expansion at 
small Q 2 entails [2TJ 




One however encounters a notable sensitivity to the perturbative behavior of 
the qq T-matrix, precisely: to the off-shell dependence of the diquark ampli- 
tudes. The cases n = 1 and n = 2 correspond to a suppression of the subleading 
amplitudes at large spacelike diquark momenta by a power of \f~P 2 and P 2 , 
respectively. Larger n implies being closer to the dominant-amplitude result. 
This sensitivity is already visible in the plain electric form factors G P E , G E , 
and it is again strong in the axial-vector diquark-photon contribution which 
depends on the diquark amplitudes via Eq. flD.27j) . In the form factor ratio, 
the systematic error band at increasing Q 2 (basically the region between full 
and dashed curve) induced by this uncertainty is quite large. Additionally, at 
the present level of sophistication the off-shell ansatz for the diquark ampli- 
tudes is just a surrogate for the correct perturbative behavior of the T-matrix 
(in terms of the bare ladder kernel which has been truncated by the diquark 
ansatz). This perturbative behavior is clearly important at large Q 2 , and even 
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if one could eventually access that region the predictive power of the find- 
ings might be questioned, at least within use of the diquark ansatz for the 
T-matrix. Such uncertainties could be removed by calculating the scattering 
matrix directly from the kernel in Eq. (J3]). 
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Fig. 17. The proton's form factor ratio for different relative behavior of the diquark 
amplitudes. Rosenbluth experimental data (circles) are taken from |113j . polariza- 
tion transfer data from [114j (squares) and |115j (diamonds). 



8 Conclusions and outlook 

The main goal of this work was to extend previous nucleon studies in a quark- 
diquark approach by implementing results obtained from Dyson- Schwinger 
and Bethe-Salpeter equations at a more fundamental level of QCD. The quark- 
diquark model can be derived from the full relativistic three-body equation 
by a series of well-defined truncations that preserve Poincare covariance. We 
have studied the nucleon mass as well as charge radii and magnetic moments 
as functions of the quark or pion mass to a point well beyond the strange- 
quark mass. Furthermore, we calculated electromagnetic nucleon form factors 
up to a photon momentum-squared of Q 2 ~ 1 — 2 GeV 2 . The former allows 
for a comparison with lattice results and their chiral extrapolations; the latter 
is of large interest in the current experimental situation regarding the ratio of 
\xGeI Gm- 

In order to link our investigation of nucleon properties to existing advanced 
meson studies we employed a rainbow-ladder truncation in the quark DSE 
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and the qq and qq scattering kernels. These correctly implement characteris- 
tic non-perturbative features of QCD such as spontaneous chiral symmetry 
breaking and yields effective confinement of the quark propagator. By adjust- 
ing the shape and quark-mass dependence of the effective coupling a{k 2 ) to 
pseudoscalar meson observables and lattice results for the quark propagator, 
meson properties are readily obtained and all parameters are fixed, yielding 
all nucleon observables as predictions of our setup. An analysis of the vector 
meson mass confirms that our results from rainbow-ladder truncated DSEs 
are within 10% of the corresponding lattice data. 

A similar situation exists for the nucleon mass, where in addition to the 
rainbow-ladder truncation a diquark ansatz is used. In principle, the diquark 
ansatz is a reasonable approximation in the infrared region where contribu- 
tions from diquark correlations are dominant. In the UV region the truncation 
causes systematic uncertainties in the electric form factors at finite Q 2 (and 
therefore also in the form factor ratio /iGe/Gm), whereas it does not affect 
the nucleon's charge radii and magnetic moments. Clearly, this problem could 
be remedied by omitting the diquark ansatz in favor of a direct solution of the 
relativistic Faddeev equation in a gluon ladder truncation, thereby advancing 
the current treatment of baryons to that of mesons in the DSE framework. 

Other issues are connected to the use of the rainbow-ladder truncation of the 
DSEs: the similar characteristics in the quark-mass dependence of both vec- 
tor meson and nucleon masses is remarakble. It suggests that effects beyond 
rainbow-ladder truncation account for the main discrepancy between the re- 
sults presented in this work and those obtained from a combination of lattice 
and chiral effective field theory methods. Part of such effects come from the 
pion cloud, i.e., the long-range qq interactions with the nucleon. The absence 
of such contributions is clearly visible in our results for the nucleon's charge 
radii and magnetic moments. 

In this respect, the interplay of pion cloud- and other contributions beyond 
rainbow-ladder truncation plays an important role. From chiral effective field 
theory and studies within the cloudy bag model it is known that pion loops are 
attractive and reduce the nucleon mass near the chiral limit, whereas its effect 
becomes small at larger quark masses. The present constant underestimation 
of the lattice data encountered in the vector-meson and nucleon masses makes 
clear that a truncation of the DSEs beyond rainbow-ladder must also provide 
further repulsive effects. 

It is conceivable that both contributions could be taken into account by di- 
rectly implementing a more sophisticated quark-quark interaction kernel in the 
relativistic Faddeev equation. Despite the need for such future improvement, 
the results presented in this work show that the quark-diquark approach is 
able to capture the major part of the nucleon's quark core, if its Green function 
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content is consistently obtained from the dynamic equations of QCD. 
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A Supplements to the quark DSE 



The quark DSE (J9j) can be reexpressed in terms of two coupled integral equa- 
tions for A(p 2 ,fi 2 ) and M(p 2 ), 

A(p 2 ,^) = Z 2 ( M 2 ,A 2 ) + S A (p 2 , / i 2 ,A 2 ), (A.l) 

M(p 2 ) A(p 2 , ^ 2 ) = M(A 2 ) Z 2 (/A A 2 ) + S M (p 2 , /i 2 , A 2 ), (A.2) 

where Ha and H m are obtained from the quark self-energy via 

S(p,/x,A) =^S A (p 2 , y u 2 ,A 2 ) + S M (p 2 , / u 2 ,A 2 ). (A.3) 

In rainbow truncation they read: 

^ / 9 9 . 9n 16nZn f d 4 q , 9 . a(A; 2 ) / 2p-kq-k\ ,. . 
M P 2 , P 2 , A 2 ) = / «M) -Li (p., + J^- J , (A.4) 

SW.P 2 ,^) = l^zlj^a.^^p-. (A.5) 

Eqs. ( 1A.HIA.2]) can be solved iteratively for chosen values of Z 2 and M(A 2 ). 
Alternatively one can employ a renormalization condition, e.g. A(fi 2 ,fi 2 ) = 1, 
and specify the current mass M (/i 2 ) at the renormalization point as an input 
value. Then both Z 2 and M(A 2 ) are determined together with A(p 2 ) and 
M(p 2 ) in the course of the iteration via 

Z 2 (fi 2 ,A 2 )= l-S A (/i 2 , / u 2 ,A 2 ), (A.6) 
M(A 2 ) Z 2 (/i 2 , A 2 ) = MGu 2 ) - E M (/i 2 , /i 2 , A 2 ). (A.7) 



Asymptotically, the DSE solution for the quark mass function reproduces the 
behavior predicted from perturbation theory: 

M(p 2 ) p2 —; ™ + ^ (A.8) 



where .FQo 2 ) = [p 2 /Aq CD J . The coefficients rh and — (qq) define the 
renormalizat ion-point independent current mass and chiral condensate. For 
finite current masses, the second term is suppressed by a factor of p 2 while in 
the chiral limit, defined by rh = 0, it determines the behavior of the asymptotic 
mass function. The renormalization-point-dependent chiral quark condensate 
is obtained from 

- (qq),, = Z 2 (fi 2 ,A 2 )Z m (^,A 2 )N c J ^ Tr D {S chiral (g, //)}, (A.9) 
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with Z m (fi 2 , A 2 ) = M(A 2 )/M(yU 2 ), evaluated at large current masses. For large 
renormalization points, the condensates are related via 

-(qq), = -(qq)H^r m - (A.10) 

For p 2 E R + , the propagator functions a v (q 2 ), cr s (q 2 ) and the coupling a{k 2 ) 
in the integrals T, A and T, M are only needed on the positive real q 2 and k 2 
axes. In this case the coupled system (1A.1UA.2]) can be solved without com- 
plications. The straightforward way to evaluate the quark propagator for a 
complex argument p 2 E C is to implement the result for q 2 E R+ and insert 
the coupling a{k 2 ) at the complex values k 2 = p 2 + q 2 — 2 p ■ q which constitute 
the region within the parabola (t±i llmi/p 2 ]) 2 , t E R+. Strictly speaking, this 
is practicable only if the integrand inside this region is free of singularities in 
k 2 , i.e., according to (1A.4I) and (1A.5I) . if the coupling is singularity-free and 
a{k 2 -> 0) -> k 4 . 
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Fig. A.l. Branch cuts in the complex q 2 domain of the quark propagator corre- 
sponding to a fixed external point p 2 and a possible integration path (see text). 

After performing all the integrations except the q 2 integral, singular points 
in k 2 lead to branch cuts in the complex q 2 plane, illustrated in Fig. IA.lt for 
instance, a pole at k 2 = generates a circular branch cut in the complex 
q 2 plane with an opening at q 2 = p 2 (dashed line). The logarithmic one- 
loop behavior (JT5l) entails that the coupling necessarily exhibits singularities 
also at k 2 ^ 0: those will generally lead to more complicated branch cut 
structures (dash-dotted line) which, however, still leave the arc q 2 = re ,arg? , 
r E R + , unharmed. A possible way to avoid all occurring branch cuts (an 
alternative method is described in [107] ) is to deform the integration contour 
q 2 E (0, A 2 ) to a complex arc that passes the point p 2 and eventually returns 
to A 2 E R in the far spacelike region. Since (J v (q 2 ) and cr s (g 2 ) must already 
be known on these complex paths, the complex DSE solution is therefore 
obtained via iteration of (1A.lHA.2j) on a family of deformed complex paths in 
p 2 . The disadvantage of this method is that successively increasing the angle 
of the respective arc is viable only up to the first pole (pair) of the propagator 
occurring in the timelike complex plane. 
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A virtue of the parameterization (Tl6|) is that these issues are largely avoided 
by its infrared behavior a(k 2 — * 0) — > k 2 , and by the large oscillations caused 
by the exponential parts which effectively shield the complex conjugated poles 
stemming from the logarithmic tail in the coupling. The remaining 1/k 2 pole 
in the Ha integral results in small numerical artifacts which are visible in the 
complex functions o~ v {p 2 ) and o~ s {p 2 ) if they are directly calculated without 
employing refined methods. However, if one employs a different coupling that 
behaves as a(0) = const, or involves more malicious singularities, these ar- 
tifacts become dominant and one inevitably has to resort to more advanced 
methods. Similar procedures can be used for evaluating meson and diquark 
amplitudes in the complex plane of the relative momentum between the con- 
tributing quarks. 



B Meson and diquark amplitudes 



B.l General decomposition 



The meson BSA has been introduced in (IPTj) together with the variables it de- 
pends on. The Dirac structure of the amplitudes is determined by the Clifford 
algebra of the Dirac 7 matrices. As for any fermion-fermion-scalar or fermion- 
fermion-vector vertex, it allows for 4 basis matrices for the pseudoscalar am- 
plitude, 

T 1 .. A (q,P) = {t,f,iif}, (B.l) 
and 12 basis elements for the vector amplitude: 

TL 12 (q,P) = r{hF,i,i?}, <f{i,r,i,if} , p*{i,r,i,i?}. (b.2) 

Implementing the negative parity requirement 



r(g,P) = - 7 4 r(Ag,AP) 7 4 , 
r"( 9 , P) = 7 4 AT"(A g , AP) 7 4 



(B.3) 



where A = diag(— 1, —1,-1,1), gives rise to the following general structure of 
pseudoscalar and vector meson amplitudes, written with full Dirac, color and 
flavor dependence: 



r(g,P) = E/r(g 2 ,^P 2 ){^(g,P)} Q/3 ® ^ w , ab , 



5 (R4) 
r ^ p ) = E/ fe v V>^^ 2 )K(g,P)} a/3 ® ® 4. 

The respective dressing functions fk{q 2 , z, P 2 ) only depend on the Lorentz 
scalars q 2 , P 2 and the angular variable z = q-P (q denotes a normalized 4- 
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vector q/^/q 1 )- By solving the meson Bethe-Salpeter equation, these dressing 
functions are obtained on the domains q 2 G R + (by use of the same methods as 
discussed in App.lAlalso for g 2 6 C), P 2 = —M 2 (i. e., on the mass shell), and 
z G ( — 1,1). Greek indices refer to the Dirac structure. The color structure 
of the meson amplitudes is diagonal, with A, B = 1, 2, 3. We are working 
with two flavors and assume isospin symmetry. The flavor matrices r\ h (a, b = 
1,2) then refer to isospin singlet (e = 0) and triplet states (e = 1,2,3); in 
the pseudoscalar case we only consider the triplet states (i.e., the pion). We 
choose the r e to be normalized to unity via Tr{r e V e } = 5 ee /, but mostly omit 
them from formulae, since they give no contribution to the Bethe-Salpeter 
equation (flTj) : for degenerate flavors the propagators in the kernel are not 
flavor-dependent. The charge-conjugated amplitudes are defined by 

t {q ,-P) = C^(- q ,-P)C-\ 
P"(g, -P) = -C r" T (-q, -P) C~ l . V ' 

where the superscripted T denotes matrix transposition. 




Fig. B.l. Notational conventions for diquark amplitudes. The quark momenta are 
q± = ±q + P/2 and the indices are Dirac/Lorentz, color and flavor indices. 

Apart from the opposite parity requirement and different color and flavor 
tensors, scalar and axial-vector diquark amplitudes exhibit the same structure 
as their pseudoscalar and vector meson counterparts. The incoming antiquark 
momentum — q_ is replaced by an outgoing quark momentum g_ which is 
reflected by the charge conjugation matrix C = 7 4 7 2 . We denote diquark 
amplitudes by the same generic symbol T: 

r(?,P) = pj s k c (q 2 ,z,P 2 )^T k (q,P)C}^ ® ^ ® s° b , (B.6) 

r%P) = Y,fr(q 2 ,z,P 2 ){iTi;( q ,P)C} af3 ® f^jp ® s^- 3 • (B.7) 

The definition of conjugation is the same as for mesons with corresponding J 
quantum numbers, see Eq. (IB. 51) . Due to the Pauli principle, diquark ampli- 
tudes must be antisymmetric under quark exchange q+ *->■ 

T(q,P) = -T T (-ql P) , (B.8) 
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where the transposition involves all Dirac, color and flavor indices. Because of 
the antisymmetry of the color anti-triplet diquark, the combination of flavor 
and spin structure must be symmetric. Therefore spin and isospin states co- 
incide for the two-flavor case: scalar diquarks correspond to an antisymmetric 
isospin singlet and axial- vector diquarks to a symmetric isospin triplet. The 
isospin singlet and triplet matrices s° b and s^" 3 are explicitly given by 

s° = (ud t - du + )/v^ = , 
s 1 = uu 1 " = (1 + <r 3 )/2 , 
s 2 = (ud f + du f )/V2 = aJy/2 , 
s 3 = ddt = (1 - a 3 )/2 , 

where o~i are the Pauli matrices, and u = (1,0), d = (0,1). The s e are nor- 
malized to unity: Tr{s e ^s e } = 5 ee i and do not give a contribution to the BSE 
integral (I23p itself. However, the color factor ~ &abe representing the diquark 
anti-triplet configuration leads to a prefactor 1/2 in front of the integral if 
the rainbow-ladder kernel is inserted: diquarks are less bound than mesons. 
The Dirac amplitudes ~ 1, 7 M in (1B.1I) and (1B.2I) are the dominant ones in a 
solution of the rainbow-ladder BSE for the lowest-mass mesons and diquarks 
and reproduce masses of the full solution within an error of < 20% [38.39.68j. 



B.2 Mesons and diquarks on the mass shell 



For the actual solution of the meson or diquark BSE, eqs. (TTTT) and ([23")) . it is 
advantageous to construct orthogonalized versions of the general Dirac basis 
elements (IB.ltlB.2j) at the respective mass pole P 2 = —M 2 . A suitable choice 
for the (pseudo-) scalar case is 

n = l, r 2 = f, r 3 = q-P$ T , T 4 = -i[f>,$], (B.9) 

where = Tp w q u and Tp V = 5^ u — P fl P u is the transverse projector with 
respect to P, and for the (axial-)vector case: 

(B.10) 

Fix - fv 

Using normalized P M simplifies the discussion of the diquark amplitudes' off- 
shell behavior (cf. App. lB.3"]) ; using normalized q^ is convenient but not further 
relevant since the diquark BSE is solved for complex q 2 in the first place. These 
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basis elements are designed such that all corresponding dressing functions 
fk{q 2 ,z, —M 2 ) are real and even in z for o = 1/2 owing to charge conjuga- 
tion symmetry. On the mass shell it is sufficient to consider 8 of the general 
12 components since on-shell vector mesons (and axial-vector diquarks) are 
transverse with respect to their total momentum P: the basis elements 9 ... 12 
in this representation correspond to the basis elements P M {1, P, $ P} which 
are purely longitudinal. Therefore, for vector mesons and analogously for axial- 
vector diquarks, one has 

r^,P)| p2 _ M2 =TJ? j:fTtf,z-M^{i7Z{ q ,P)\ } . (B.ll) 

vc k = l a P 

The orthogonality relations for the basis elements are 



where 
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Tjf Tr{7f rj] 
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T = - 63 = ? 
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(B.12) 



l-z 2 , (B.13) 



Contraction of ifTTj) and (123D with the basis elements and using the orthogo- 
nality relations leads to coupled homogeneous integral equations for the co- 
efficients fk(q 2 , z, —M 2 ). Since they depend only on Lorentz-invariant scalar 
products, these equations can be solved in any arbitrary frame. Exploiting the 
0(4) symmetry of the problem, one can employ a Chebyshev decomposition 
for numerical convenience: 

n 1 

f k (q 2 } z,-M 2 ) » fk(Q 2 )U n (z), / dzVT^U m (z)U n (z) = - 5 mn , 

n=0 I 1 

(B.14) 

where the U n (z) are Chebyshev polynomials of the second kind. For a suitable 
choice of basis elements, only a few Chebyshev moments fk{q 2 ) have to be 
taken into account to match the full solution [3"8~|I39"] . We will use this observa- 
tion below and reduce the angular dependence to a constant, see App.|E]and 
Sec. [71 By introducing an artificial parameter A(P 2 ) at the right-hand side of 
the homogeneous BSE (123]) . the equation becomes an eigenvalue problem for 
A(P 2 ) where a bound-state solution is obtained for A(— M 2 ) = 1. 

As mentioned in Sec.[2j bound-state amplitudes can be normalized via a 
"canonical" normalization condition, Eq. (151) , derived from the pole condition 
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([2]). In the case of a two-body Bethe-Salpeter amplitude obtained via rainbow- 
ladder truncation, the ladder kernel is independent of the total momentum P 
such that only the derivatives of the propagators contribute. For pseudoscalar 
and vector mesons this reads: 
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dP 2 



Q(P 2 
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dP 2 



Qt{P 2 



(B.15) 



with 



Q itMu) (P) ■= I ^Tr D {rM( q ,-K)S( q+ )rM(q,K)S(-q„)} 

J \ Z7[ ) K 2 =-M 2 

(B.16) 

The notation Q^ u ^ covers either the pseudoscalar quantity Q or its vector 
counterpart Q^, where Qt = Tp V Q^ u /3 is the transverse component of the 
latter. Tr^ denotes a Dirac trace; the color-flavor trace is 1 since we use nor- 
malized color and flavor matrices. The Bethe-Salpeter equation ( fI7|) . together 
with the normalization ( lB.15llB~T6l) . completely determines the meson ampli- 
tudes on the mass shell. The pion decay constant is calculated via 
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-M 2 



(B.17) 

where the prefactor J Nc / Np is again a consequence of the color-flavor nor- 
malization. We note that Np = 2 here, since we are dealing with bound states 
of two light (anti)quarks here. 



The normalization condition for the diquark case can be obtained by differ- 
entiation of Eq. (124"]) with respect to P 2 at the mass pole, insertion of the 
diquark BSE (K^ 1 — Go) T = and using the fact that the ladder kernel K 
is independent of P 2 (G = —^S S contains a symmetrization factor 1/2 for 
the diquark case 
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Qt(P 2 ) = i 



P 2 =-M 2 



<B.m 



with 



Q ( " V \P) ■= \ J ^Tr D {T^\q,-K)S(q + )T^(q,K)S T (q_)} 



(2k) 



K 2 =-M 2 

(B.19) 
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B.3 Off shell ansatz for the diquark amplitudes 



With the methods outlined in this work we cannot gather any information on 
diquark amplitudes for off-shell total diquark momenta P 2 ^ —M 2 . On their 
mass shells, the amplitudes are obtained from the Bethe-Salpeter equation 
(1231) and equivalent to the residues of a scalar or axial-vector quark-diquark 
vertex at the lowest-lying diquark poles. An inhomogeneous diquark BSE, to 
be solved in an analogous way as an inhomogeneous meson BSE [23], also pro- 
vides only very limited information on the off-shell behavior of these residues. 
Therefore we choose a reasonable ansatz for the diquark's off-shell behavior 
based on general assumptions. 

We start from Eqs. (lB.6tlB.7l) and for convenience discuss the offshell depen- 
dence of the diquark amplitudes in terms of the basis elements Tk{q,P) and 
T%(q,P) while the dressing functions are left unchanged at their mass-shell 



values: f k (q 2 ,z,P 2 ) = fl (q 2 ,z,—M 2 c / a:v )). Usage of the orthogonal on- 



shell bases ( 1B.9I) and ( IB.lOj) . where each appearance of P^ has been nor- 
malized, entails that the correct behavior at the mass shell and for P 2 = 



therein, or generally by an arbitrary function h(P 2 /M 2 ) satisfying h(—l) = 1 
and h(0) = 0. This then also applies to each transversal projector: Tp V — > 
5^ u — h 2 (P 2 /M 2 ) Vff . Our ansatz for the full Dirac part of the diquark ampli- 
tudes therefore reads 



can be guaranteed by attaching a factor 




i 



T sc (q,P) 



Y,f S k C (q 2 ,z,-Ml) ll 5 T k (q,P)C 



(B.20) 



k=l 
8 



Y.fu{Q 2 ^-Ml)ir^P)C 



(B.21) 



k=l 



where the scalar basis is given by 



n = g sc i 




(B.22) 




r 4 = -i h sc [f, jf] , 



and the axial-vector basis by: 




(B.23) 
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In addition to the functions /i sc (av) : = M x sc( av )) = h(P 2 /M^.,^), which so 
far are arbitrary except for the two above requirements, we also attached a 
function g sc ( av ) := g(x: sc {a-v)) with g(—l) = 1 to each of the two leading ampli- 
tudes in order to be able to modulate them separately from the others. The 
transversality condition on the mass shell for the axial- vector amplitudes need 
not be stated explicitly since it is already ensured by the transversal pole in 
the axial-vector diquark propagator: each diquark amplitude in the subsequent 
calculations appears in conjunction with the respective propagator. Likewise, 
the purely longitudinal components related to 12 (g, P) off the mass shell are 
generated by the longitudinal contribution to the diquark propagator which 
is suppressed by a factor of P 2 + M 2 V on the mass shell. 

The asymptotic behavior for P 2 — > oo of the diquark amplitudes is provided 
by the limits g(x — ► oo) and h(x — > oo). We construct these functions based 
on two assumptions: firstly, all subleading amplitudes related to Tfc(g, P) and 
r£(q, P), k > 1, should be suppressed by a power of at least \f~P 2 compared 
to the leading ones Ti(q,P) and rf(g, P) such that the diquark amplitude's 
perturbative limits are indeed 27 s C and ry M C. Secondly, the dominant ampli- 
tudes should behave like y/P 2 in the ultraviolet to ensure via Eq. (124"]) that the 
diquark propagators asymptotically behave as D(P 2 — > 00) — > 1/P 2 . Simple 
parameterizations that fulfill all of these requirements are 

g(x) = Vx~+2, h{x) = ~^^^, n>\. (B.24) 

For large n the subleading amplitudes are suppressed at spacelike momenta 
and provide support only in the neighborhood of the mass shell. This resembles 
the case where only the dominant diquark amplitudes have been taken into 
account. Furthermore, a global function depending on P 2 which is attached 
to all scalar or axial- vector amplitudes does not change the product T DT 
since it would also appear in the diquark propagator and leave the T- matrix 
itself (and therefore baryonic observables) invariant. For instance, one could 
divide all diquark amplitudes by g(x) such that the leading amplitudes become 
constant in the ultraviolet while all the others are asymptotically suppressed: 
then also the diquark propagator would become constant for P 2 — > 00. This 
has been applied, e.g., in Ref. [20J. 



B.4 Diquark propagator 

The diquark propagator is obtained by Eq. (|2~4"j) for all values of the squared 
diquark momentum P 2 . This equation can be written as 

(£)-!) (/^)(p) = K^ U \P) + Q^ V \P) , (B.25) 
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where K^ v > represents the first term involving the inverse kernel and Q^ u ^ 
the quark loop term (in contrast to Eqs. (IB.16I) and (IB.19I) . now without the 
P 2 dependence of the diquark amplitudes held fixed), 

Q iMu) (P) ■■=\j-^^ D {f^( q ,-P)S( q+ )T^( q ,P)S T ( q ^} . (B.26) 
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Fig. B.2. Graphical representation of the diquark propagator of Eq. (|B.25D , 

K^ V \P) involves a 2-loop integral, to be evaluated for complex momenta. We 
circumvent its explicit calculation by observing that the pole and unit residue 
conditions (in the axial-vector case only for the transverse part) 



-M z 







-M z 



[B.27) 



where the derivative ' denotes d/dP 2 and each quantity F^ u = (D 1 ) /w/ , K^ v 
and is decomposed as 

F^(P) = F T (P 2 ) Tp U + F L (P 2 ) L P V , (B.28) 

with F T = TjfF^/3, F L = L P V F» V , L P V = P»P U and Tjf = 8^-L P u , impose 
the following constraints on K(P 2 ) and Kt{P 2 ): 

K {T) {-M 2 ) = -Q (T) (-M 2 ) , (B.29) 
K[ T) (-M 2 ) = -Q{ T) (-M 2 ) + 1 . (B.30) 

Asymptotically (i. e., for P 2 — > oo) K ^ becomes dominant in (1B.25|) since the 
ladder kernel is independent of the diquark momentum whereas Q^ u ^ contains 
the product of two quark propagators which behaves as 1/P 2 . The ansatz 
(!B~24l entails K ^ -> P 2 and Q^ u) -> const. Analyticity at P 2 = requires 
Kt(0) = Kl(0)] the corresponding relation Qt(0) = Ql(0) is guaranteed 
by the definition ( IB. 261) . We employ a parameterization based on the above 
restrictions: 



I 2 ) 

L scJ ) 



K(P 2 ) = (l - Q'(-M 2 C )) (P 2 + M 2 C ) - Q(-Ml f ) oi 

K^(P) = {(l - Q'A-M 2 J) (P 2 + M 2 av ) - QA-Ml)} ST . 

With Eqs. flBT20]) — (|BT26|l and Eq. flB~31~j) . the off-shell behavior of the T-matrix 
is completely determined. 
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C Quark-diquark amplitudes and quark-diquark BSE 



The matrix- valued Bethe-Salpeter amplitudes (p, P) which were intro- 
duced in Sec. [5] feature a decomposition constructed from the same Dirac basis 
elements as used in the meson and diquark case: 



k=l 1 V<3 



(C.l) 



tab 



Here p is the relative momentum between quark and diquark and P is the 
total nucleon momentum on the mass shell: P 2 = — Mjy. The positive parity 
condition for the full Faddeev amplitude translates into positive parity of the 
quark-diquark amplitudes. The constraint of positive energy for the nucleon, 
expressed by the positive-energy projector 

A+OP) = \U+ -TF-) , (C2) 



2 V iM 



N 



halves the number of possible Dirac basis elements via Jf A + = iM^A + from 
4+12 to 2 Dirac matrices {1,^} for the scalar quark-diquark amplitude and 
6 matrices {7^, 7^, p^p 1 , P M , P^p 1 } for the axial-vector one. A possible or- 
thogonal basis set is T\ — 1, Ti — —ifr for the scalar and 

rf = 7^ t% = % (7^ j5 T - p£) 

T2 1 = ^ r 5 M =PtPt- y7r ( C -3) 

r 3 M = P M r 6 " = zP^t 

for the axial- vector basis. Similarly as before, pj. = Tp U p u , where Tp V is now 
the transverse projector with respect to the total nucleon momentum. The 
orthogonality relations are 

^Tr{r,r i A + } = S ija f(z), -Tr^/V ^7 5 A + } = S^z) (C.4) 



and 
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2(1 -z 2 ), (C.5) 
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Particular linear combinations of the eight Dirac structures Tfe(p, P) A+(P) 
and T%(p, P) 7 5 A + (P) lead to a partial wave decomposition in terms of quark- 
diquark total spin and orbital angular momentum eigenstates in the nucleon's 
rest frame [rBfllb'j . These combinations correspond to relative s, p and d waves 
which implicates that the nucleon is not spherically symmetric. 

The flavor matrices t e in the quark-diquark amplitudes, when applied to proton 
and neutron isospinors (1,0) and (0,1), yield the "quark remainders" of the 
full proton and neutron flavor wave functions constructed from the Clebsch- 
Gordan prescription after removing the diquark contributions. They are given 
by t° = 1, t 1 = (<ti - ia 2 )/VQ, t 2 = -a 3 /V3 and t 3 = -{a x + i(T 2 )/VE- 
The flavor traces of the quark-diquark Bethe-Salpeter equation for proton 
and neutron are obtained by projection on the right and the left with (1,0) 
and (0, 1), respectively. Similarly, the antisymmetric color factor of the full 
Faddeev amplitude is obtained as product of diquark and quark-diquark color 
factors and has its origin in the diquark amplitudes ( 1B.6UB.7l) . 

The quark-diquark amplitudes are normalized by an analogous normalization 
integral as given in Eq. ([5]): 

5 -(^. = 1 . (c . 6) 

In this context, $ is the quark-diquark amplitude as solution of the Bethe- 
Salpeter equation (I26I) . K is the quark-diquark kernel (127]) and Go the product 
of dressed quark and diquark propagators. The normalization condition is 
equivalent to the normalization of the electric charge of the proton, G P E (0) = 1 
[T7] . In contrast to the meson and diquark case, Eqs. (1B.15P and (1B.18I) . the 
quark-diquark kernel depends on the total nucleon momentum P such that 
dK~ l /dP 2 cannot be omitted. 



D Construction of the electromagnetic current diagrams 



The electromagnetic current of Eq. (1331) has the explicit form 

(D.l) 

where Pj and Pf = Pi + Q are incoming and outgoing on-shell nucleon mo- 
menta. The loop momenta Pi and pf are arbitrary, a, (3 = 1 . . .4 are quark 
and a,b = 1 ... 5 are diquark indices. The quark-diquark amplitudes $ a are 
the solutions of the quark-diquark Bethe-Salpeter equation ( 1261) . The quantity 
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X M ' is given by 

X », ab = X p* (27r) 4 J4 {pf _ pt _ {1 _ 7])Q) + 

+ X^ b (2tt) 4 5 4 ( Pf -p t + V Q) + Xr b , 1 ' J 

with 

= {s( P+ ) i^(p + ,p-) 5-(p_)} Q ^° 6 (p d _) , (d.3) 

(x dq )^f = ^(p_) {i^y rgfW.i*-) D b '\ Pd j)} , (d.4) 

(X K )^ = l?«'(p^){5(^)Jif^(p /) p ( ,P / ,P 4 )5(p_)}^ P 6 'W) . 

(D.5) 




Fig. D.l. The five diagrams that contribute to the nucleon's electromagnetic cur- 
rent, corresponding to Eqs. (|D.3|) (|D.9p . 



The first two contributions correspond to the impulse approximation and con- 
tain quark-photon and diquark-photon vertices. The third term is a two-loop 
diagram and reflects the gauged quark-diquark kernel. The diquark propagator 
and therefore also the quark-photon coupling X^ ab is either scalar (a, b = 5) or 
axial- vector (a,b = 1 ... 4) while X^ and Xjj? a mix scalar and axial- vector 
quantities. The quark and diquark momenta are: 

P-=Pi + i]Pi p d _ = -pi + {l-rj)Pi , 

P+=Pf + V p f Pd+ = ~Pf + (1 - v) p f ■ 



The gauged kernel contains the exchange quark diagram and the seagull con- 
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tributions: 



with 



= r^(p!,p d _) g) 5(g)} ffo, -p d+ ) (D.7) 

K$g> = M^ b (k uPd _, Q) S T (q') T a (p 2 , -p d+ ) (D.8) 
K^ b = r b dq ( Pl ,p d _) S T (q) M^ a (k 2 , -p d+ , Q) (D.9) 



and momenta: 



p+-q P-~q' 
q = p d --p+ , pi = — - — , p 2 = — g 

1 — Pd+ — V- > Kl = o ' fc 2 = = ' 



For explicit calculations we work in the Breit frame where 



^t^ = {o, 0, 0, iMVTT^} , Q = {0, 0, |Q|,0} , (D.10) 



with r = Q 2 /(4M 2 ). The electromagnetic current is completely determined 
by specifying the quark-photon vertex T^, the diquark-photon vertex r^ q , and 
the diquark-quark-photon vertex or seagull M M . 



D.l Quark-photon and diquark-photon vertex 



With the notation P = {p+q)/2 and Q = p—q, where p and q are outgoing and 
incoming quark or diquark momenta, the general form of the quark-photon 
vertex and the scalar and axial-vector diquark-photon vertices is 



12 



I>, q) = E /*V, ^ Q 2 ) (r q )£ (P, Q) , (D.11) 
fc=i 

r^, ?) = E /r d V, ? 2 , Q 2 ) (P, 0) , (D.12) 
k=i 

rgfV ?) = E tfTV. 9 2 , Q 2 ) (r dq r k a " (P, Q) , (D-13) 
fc=l 
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with the basis elements: 



fa)Li...u ( P > Q) = 7 M {l,$,f,f0}, (D.14) 

P"{1,#, P,P 

Q"{i,0, p, p m 

(WLT.iJ^) = n^^.nQ ,^}, (D.i6) 

P^{<5 Q/3 , P a pP, Q a Q t3 , P a Q p , Q a P 13 } , 
Q^{5 aP , P a P 13 , Q a Q p , P a Q p , Q a pP} . 

All vertices satisfy vector Ward-Takahashi identities which reflect electromag- 
netic current conservation. They constrain the longitudinal contributions ~ 
by relating them to the corresponding quark or diquark propagators: 

Q*T*(p,q) = S- 1 (p)-S-\q) (D.17) 
Q^T^p, q ) = D-\p 2 )-D-\q 2 ) (D.18) 
Q" r£?V q) = D-}{ P ) - D~l{q) (D.19) 

Implementation of the WTI for the quark-photon vertex with the general 
quark propagator (|SJ) leads to the expression 



8 (D.20) 
+ AB(p 2 , q 2 ) (p + g)« + lg" £ / fe V, g 2 , Q 2 ) (r^ (P, Q) , 

fc=i 

where B(p 2 ) = M(p 2 )A(p 2 ) and 

i.e., a sum of the Ball-Chiu vertex |117j and 8 unknown transverse compo- 
nents. The scalar diquark-photon vertex then reads 

K q (P, q) = ( P + qY A2T V, q 2 ) + /rV, q\ Q 2 ) Tq P u (D.21) 
and the axial- vector diquark-photon vertex: 

Kf(p,q) = (p + qY{AD T \p 2 ,q 2 )6^-Aa D (p 2 ,q 2 )p a qP} 

- {v D (p 2 )^p a + VD(q 2 )^ a q^} (D 22) 

- Tq e /r d v, q 2 , q 2 ) (ra^r" & q) , 
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with <Jd(p 2 ) '■= [D^iP 2 ) — D l 1 {p 2 )^j jp 2 such that the inverse axial-vector 
diquark propagator is written as 

D~p{p) = D T \p 2 ) 5 a ? - o D (p 2 ) p«pP • (D.23) 



The transverse parts of the vertices (1D.20l4D.22l) must vanish for Q 2 — > 0, 
either by a Q^ L dependence of the basis elements or by a vanishing amplitude 
at Q 2 = 0. This guarantees the differential Ward identities 



r£(p,p) 

r dq(P,P) 



dp^ 
dD- l (p 2 ) _ 
dp^ 

dDapiP) P 2 



iYA(p 2 ) + 2p^ if 



. . dA dB\ 



dp 2 ^ dp 2 I 



2 2 



2p», 



d 



dp^ 



dp^ v 



(D.24) 
(D.25) 
(D.26) 



In the axial- vector case the symmetry requirement T^^(p,q) = T^ a (q,p) 

entails f 3 = -f u f A = f 2 , fg = -f%, and /i, 2 ,5,6,7,8 = /i,2,5,6,7,8 where 
fi(p 2 ,q 2 ,Q 2 ) = fi(q 2 ,P 2 ,Q 2 ), therefore 6 unknown dressing functions are left 

Eng. 






m)= + =(W 





Fig. D.2. The full diquark-photon vertex. The first two rows correspond to 
Eq. (|D.27j) . The last row can be dropped for the scalar-axialvector transition ver- 
tex; for scalar and axial- vector diquark-photon vertices an ansatz is employed, see 



In the actual calculation we follow Ref. [20] and express the diquark-photon 
vertices by the "gauged" diquark propagators, cf. Eqs. (1B.2511B.26l) and the 
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discussion in Sec.|6j 



fd q a6 (p,g) 



d A k 



Tr^ ab (p,q,k)+ (D.27) 
+ T a (k, -p) S(k' + ) M^(k, q, Q) S T (ki) + 
+ M^ a (k, -p, Q) S(k+) T b (k, q) S T (k_) 

with the momenta k± = ±k + q/2 and k'± = ±k + p/2. The quantity 

Af ab (p, q, k) =T a (k- Q/2, -p) S(k+) T b (k, q) x 
S(k-+Q) T»(k„+Q,k„)S(k_ 



(D.28) 



denotes the quark loop where the photon couples to the upper quark line (see 
Fig.[EL2]). The seagull contributions M^ b and M fta that appear in Eq. (11X271) 
are defined in App. ID.21 They are necessary to guarantee the Ward-Takahashi 
identities for the scalar (a, 6 = 5) and axial- vector (a, b = 1 ... 4) vertices 
and transversality for the scalar-axialvector transition vertices: Q 11 ^dq b = 
Q^T^ 5 = 0. To show this one uses the relations T T (k,P) = —T(—k,P) } 
T^ T {k,P) = r M (— k,P) for the Dirac parts of the diquark amplitudes which 
follow from the general antisymmetry relation (IB. 81) . Eq. ( 1D.27I) furthermore 
ensures the correct symmetry behavior T^ b (p, q) = —T^ b5 (—q, —p). The color 
and flavor factors in Eq. (1D.27|) have already been worked out. They entail 
e_ = e + = ed q /2 = 1/2 in the scalar and axial- vector case and e_ = — e + = 
1/2, ed q = for the transition vertex (e± and appear in the definition of the 
seagulls, cf. App lD.2l ). The color trace for all contributions is 1. For instance, 
using the diquark flavor matrices Si (i — ... 3) of App. IB.ll entails that the 
quark-loop contributions of Eq. (ID. 271) are given by 



Tr {sf Sj Q} Af ab + Tr {sf Q Sj } Af ab = 2 Tr {s^ Sj Q} Af 



ab 



(D.29) 



where Q = diag(g n , qd) is the quark electric charge matrix, and Aj_ denotes the 
counterpart to Eq. (1D.28|) where the photon couples to the lower quark line. 
In App. lD.3l the diquark charge factors 



2Tr{sf Sj Q} 



V 



q u + qd 





Qd - q u 


\ 





2q u 








Qd - Qu 





Qu + Qd 














2gd / 



(D.30) 



are explicitly attached to the current matrix diagrams at each occurrence of 
the diquark-photon vertex. 



In order to satisfy the Ward-Takahashi identities ( lD.181lD~T9l) for the full di- 
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quark propagators (1B.25j) . also the diagram involving the gauged inverse ladder 



kernel has to be taken into account when constructing the scalar and axial- 
vector diquark-photon vertices. It depends on the diquark momentum only via 
the diquark amplitudes, therefore its generic form is r a K~ 1 M fJi,b +M fl,a K^ 1 ^. 
Since we employed a parameterization for K^ u ' we have to choose an ansatz 
for this contribution. We use the generic vertices of Eqs. (ID.21HD."22l) . drop the 
transverse parts and replace D^ 1 — > K, ap —>■ ctr--i. With the ansatz (IB. 311) 
for K the full diquark-photon vertices read 

q) = f q) + (p + qY AK(p 2 , q 2 ) , 
r!Sfip, q) = fg*(p, q) + (p + qT AKrip 2 , q 2 ) ST* . 



D.2 Seagulls 



(D.32) 



Seagull contributions represent the photon coupling to the diquark amplitudes 
and reflect the diquark's internal substructure. The Ward-Takahashi identities 
for the seagulls are given by [17)119] : 

Q»M»(q,P,Q) = +e.{T(q + ,P)-T(q,P)} 
+ e + {T(q^P)-T(q,P)} 
-e dq {T(q,P + )-T(q,P)} , 
Q"M»(q, -P, Q) = - e + {f{q + , -P) - f (q, -P)} 
-e_{r(g_,-P)-f( 9 , -P)} 
+ e dq {T(q,-P_)-T(q, -P)} , 

where q is the relative momentum, P is the diquark's total momentum, Q is 
the photon momentum, q± = q ± Q/2, and P± = P ± Q. The charges e+, e_ 
and edq correspond to quark and diquark legs (cf. Fig. ID.3j) . For Q — > the 
WTIs reduce to the differential Ward identities: 

MW) = (^^- e ,jLy {q , P) , ^ 

M-( 9 ,-P,0)^(^±A_ e M f(9 ,-P). 



To simplify the discussion, we rewrite the diquark amplitudes (IB .2 014131231 in 
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Fig. D.3. Conventions for the seagulls. In accordance with the Ward-Takahashi 
identity (1ET321) . out going charges are taken to be positive, incoming charges as 
negative. 

terms of the general basis elements r]f (q, P) of Eqs. (1B.1HB.2"]) . i. e. 



r sc (g,P) = £ (g^, P 2 ) i T 5 <r fe (g, P) C 



k=l 



k=l 



sc 
k 



where z = q ■ P. The scalar and axial-vector diquark dressing functions f l 
and f k v are the respective linear combinations of the coefficients in Eqs. ( lB.20f - 
IB.21I) and therefore carry now a dependence on the total diquark momentum 
P 2 . That being said, we also drop the superscripts 'sc' and 'av' for the pur- 
pose of brevity. A possible ansatz for the seagulls corresponding to the scalar 
diquark amplitudes that satisfies the WTI and is free of singularities is then 
given by 

-z 7 5 M^(q,P,Q)C^ = 

= E {e-(V + T k r k (q + , P) + e + {V^ k r*(g_, P) - e dg V» r k (q, P + )} (D . 34) 

+ ^pV(/3 + Uf) - e dq (f2 + fdh" , 
and the corresponding ansatz for the axial-vector seagull contribution is 
-iM^' a (q,P,Q)C^ = 

= £ \e.(V + r k r k (q + , P) + e + (V_)£ r«(g_, P) - e dq V» r«(q, P + ) 
k=i 1 

+ ^{7Y(/3 + Uf) + ^ a (f 5 + hf) + {5^i + q a Y)(f7 + fsf)} 
e_ + e+ 



+ 



QVC/V + fsf) ~ e dq \ l a (f2 + hi) + q a (h + hi) \Y 



(D.35) 
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The conjugated seagulls are obtained from 

M»(q, -P, Q) = -C(M») T (-q, ~P, Q) C- 1 , 
M^ a (q, -P, Q) = C(M^ a ) T (-q, -P, Q) C' 1 . 

The quantities (V±)£ and VJf are defined by 

(V ± r k = X£ {f k (ql z, P 2 ) - f k {q\ z, P 2 )} 

V^=^{/ fc (g 2 ,z,P 2 )-/ fe (g 2 ,^P 2 )} 

q»-2aq-PX» r 2 „ 2 2 2 



rM _ (g±Q/4)" - = {P + Q/2Y 
± {q±Q/A)-Q> ~(P + Q/2)-Q' 



z±=?±-P, ^ = g-P+, l±<7± 



1 + (T = Vp^- 



In the case where only the dominant diquark amplitude is retained and its 
dressing function f\ (either scalar or axial- vector) depends solely on q 2 , these 
ansatze reduce to the forms M M = W^i^C and M^ a = W^i^C (c.f. 
Refs. [17115120121] ). where 

W» = e„X» (h{ql) - h{q 2 )) + e + X» - A(g 2 )) . 

The seagulls vanish if the diquark amplitudes are taken to be pointlike, i.e., 
fi(q 2 ) = const. 



D.3 Color, flavor and charge coefficients 



The current matrix diagrams still have to be endued with color and flavor- 
charge coefficients. The color traces for the impulse approximation and ex- 
change/seagull diagrams are given by 

$ba ^ab _ , &ba £aed £ceb £>cd _ -, / n o K \ 

73 73 ' 7s 7/1" 7/T 7s ' { } 

respectively. With the diquark and quark-diquark isospin matrices of Secs. IB.il 
and Si and ti, i = ... 3, and the quark charge matrix Q = diag(g u , qj), 
the flavor-charge matrices for the quark-photon, diquark photon and exchange 
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diagrams read: 

E^tjQtj, £t?tj2Tr{sf Sj Q} , £t? Sj Qs?tj. (D.37) 



The traces for proton and neutron are obtained by sandwiching these expres- 
sions between the isospinors u = (1,0) or d = (0,1), respectively. The index 
range of the sums in (ID. 371) depends on the type of the quark-diquark ampli- 
tude in the initial and final state: e. g., for a scalar quark-diquark amplitude in 
the final state (i. e., on the left-hand side): % = 0, for an axial- vector amplitude: 
i = 1 ... 3. For instance, with the four contributions: 



2 



S^S: ut(tjs QsSto)u 
S«-A: utfEtJsiQ4tj)u = -^^, 



S 



U t ^.ifsoQsft^u 



3 3 



A-A: utl^EtfsjQsftjju 

, i=l j=l 



2q u + qd 
2v/3 

4g M - g c i 
6 



the full exchange contribution to the proton's electromagnetic current (includ- 
ing the color factor) is given by 

= -f J% + 2j ^f {J% + J%) + JTa , (D.38) 

where Jfjg, Jf* , Jf^ and Jf* denote the Dirac parts of the exchange contri- 
butions in Eq. (ID.lj) . For the seagull contributions, the coupling of the photon 



to both quark lines and the diquark line in the seagull amplitudes (cf. Fig. ID. 31) 
have to be taken into account. Therefore all occurences of e_, e + and Cd q have 
to be replaced by the combined flavor-charge factors 

SG: e.-^tfsjQsftj, e+^^ttQsjsttj, 

ij ij 
ij 



e ^ ^ 51 ^ s i s i 2 Tr { s i s i Q} 

ij 



and equipped with the overall color factor — 1. 
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E Singularities and resulting boundaries 

Each of the momenta (1251) that enter the quark- diquark BSE can generically 
be written as 



with loop momenta qi (q 2 G R+) and the nucleoli's total momentum P 
(P 2 = —M 2 ). The real coefficients a x , (3x depend on the momentum par- 
titioning parameter r]. The complex argument l 2 x of each quantity X (quark 
propagator, diquark propagator and amplitudes) that is sampled by the BSE 
is therefore bounded by the parabola (t ± i(5xij]) M) 2 , with t G R + . Since we 
do not explicitly include residue contributions in our calculations, the largest 
possible parabola is the one which contains the singularity (if any) in the 

quantity X closest to the origin at l\ = p x : (t±im x ) 2 , with m x = |I m \/pi^|- 
For instance, the quark propagator could have complex conjugate poles, the 
diquark propagators by construction exhibit timelike poles (e. g., for the scalar 
diquark propagator: p 2 x = —m 2 c ) etc. For each momentum l x , this leads to a 
restriction M < mx/\f3x(v)\) an d in total to an upper limit for the nucleon 
mass: 



Here m q is the quark "pole" mass (cf. Table [2]), = min(m sc , m av ) is the 
smallest diquark mass (i. e., the scalar diquark mass) and A is the "pole" mass 
of the closest singularity in the relative momentum entering the diquark am- 
plitudes. The dependence of the diquark amplitudes on the total diquark mo- 
mentum is modeled via Eq. (1B.24j) to exhibit its first singularity at the diquark 
poles and therefore identical to the condition involving m^. If any of the above 
quantities is free of singularities, its contribution can be removed from the 
bracket in (1E.2[) . As mentioned in Sec.EJ the momentum partitioning parame- 
ter r) can be chosen to maximize the upper limit f{rj) for the nucleon mass. The 
analogous analysis for the quark momenta in the meson and diquark BSEs with 
momentum partitioning parameter a entails f(a) = min { m q /<j, m q / (1 — a) } 
with the maximum f(a = 1/2) = 2m q . This justifies the choice o = 1/2 and 
leads to the restrictions m mes0 n, ^diquark < 2m g . 

With qi = Q and the replacement P 2 = —{M 2 + Q 2 /4), Eq. (1E.1I) can also be 
applied to the momenta appearing in the form factor diagrams. If the nucleon 
mass is already fixed, this leads to an upper limit for the potential photon 
momenta: Q 2 < 4(/(r/) 2 — M 2 ), which is typically of the size of 1 — 2 GeV 2 and 
explains the restrictions encountered in Sec. [71 Possible additional singularities, 
e.g., in the inverse quark propagator's dressing functions A and M (which 
enter the Ball-Chiu vertex) have to be considered in this case; those may 
occur in a parametrization where a v and a s are entire functions [2"CT|21lll20j . 




(E.l) 




(E.2) 
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One can choose r\ and f{rj) for each current diagram in Fig. ID.il separately 
if for each value of rj a corresponding quark-diquark BSE is solved, since the 
quark-diquark amplitudes depend on the momentum partitioning parameter 
but the single current contributions must be independent of 77. 

The diquark amplitudes of Eqs. (lB.6l - IB~7l) depend on three complex variables 
q 2 , z, and P 2 , where possible singularities in q 2 and P 2 are avoided by virtue of 
Eq. flE.2j) . The diquark BSEs fl23|) are explicitly solved for z G (—1, 1) in terms 
of a Chebyshev expansion, Eq. flB.14j) . Within the convergence region \z\ < 1 
of the Chebyshev polynomials this solution can be analytically continued into 
the complex z plane. The quark-diquark BSE however samples the diquark 
amplitudes also in the region \z\ > 1 which is not accessible via a Chebyshev 
expansion. A BSE solution for all z G C would require refined methods due 
to the intricate singularity structure for \z\ > 1; furthermore, the residues of 
those singularities would have to be taken into account when inserting the 
diquark amplitudes into the quark-diquark BSE. We circumvent this problem 
by retaining only the zeroth Chebyshev moment in each diquark amplitude, 
i.e. assuming no z dependence at all: fk(q 2 , z) ~ fk(q 2 ) in Eq. flE.2j) . This is a 
reasonable approximation at least for \z\ < 1 where the full solution is known 
due to the strong suppression of higher Chebyshev moments. 
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